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A deep narrow-band survey for Lya emission carried out on the VLT-FORS has 
' revealed 98 Lya candidates down to a flux limit of FLyQ ^ 4 x 10^^^ erg s^^ cm^^ 

r \ , in a volume of 5500 comoving Mpc'^ aX z = 2.4 centered on the hyperluminous quasar 

' HE0109-3518. The properties of the detected sources in terms of their i) equivalent 

' width distribution, ii) luminosity function, and iii) the average luminosity versus pro- 

Qh, jected distance from the quasar, all suggest that a large fraction of these objects have 

Q ■ been fiuorescently "iUuminated" by HE0109-3518. This conclusion is supported by 

\ comparison with detailed radiative transfer simulations of the effects of the quasar 

. illumination. We therefore have a unique opportunity to directly detect and image in 

Ci ' emission dense gas, independent of any associated star formation activity. 18 objects, 

a much larger fraction than in "blank-field" Lya surveys at similar redshifts, have 
a rest-frame Equivalent Width (EWq) larger than 240 A, the expected limit for Lya 
emission powered by Population II star formation. 12 sources among these do not 
, have any continuum counterpart in a deep V-band imaging of the same field (reaching 

l/^ ' to V{la) ~ 30.3 AB). For these, a stacking analysis indicates EWq > 800A (Icr), 

, effectively ruling out Lya powered by internal star formation. These sources are thus 

•/^ ■ the best candidates so far for proto- galactic clouds or "dark" galaxies at high-redshift, 

\ whose existence has recently been suggested by several theoretical studies. Assuming 

, they are mostly ionized by the quasar radiation, we estimate that their gas masses 

^SJ ■ would be about Mgas ~ lO^M© implying that their star formation efficiencies (SFEs) 

are less than 10~^^ yr"^, equivalent to a gas consumption time of 100 Gyr, several 
times below the SEE of the most gas-rich dwarf galaxies locally, and two hundred 
times lower than typical massive star-forming galaxies at z ^ 2. We have also dis- 
^ , covered extended, filamentary gas, also likely illuminated by the quasar, around some 

H ■ of the brightest continuum-detected sources with EWq > 240A. The morphology and 

luminosity of this extended emission are compatible with the expectations for cir- 
cumgalactic cold streams, as predicted by recent numerical models, but other origins, 
including tidal stripping, are also possible. 
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1 INTRODUCTION 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile un- 
der program 085.A-0171A. 
t E-mail:cantal@ucolick.org 



The Intergalactic Medium (IGM) plays a crucial role in the 
formation and evolution of the structures in the Universe. 
Galaxies form from the IGM that settles into gravitational 
potential wells deep enough for the gas to collapse (Rees & 
Ostriker 1977, White & Rees 1978). Radiative cooling dis- 
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sipates thermal energy allowing the gas to concentrate at 
the center of the halo and to generally form a (dense) disk. 
Following fragmentation, star formation takes place. Most 
of the observational studies conducted so far have captured 
the galaxy formation process from this stage on. Despite 
their importance, the early phases which may involve po- 
tential wells too shallow to efficiently form stars are largely 
unobserved at present and remain poorly constrained. Key 
questions about how the IGM is converted into stars at high- 
redshift and how efficient this process is remain unanswered. 

Recent theoretical studies have suggested that gas-rich, 
low-mass haloes (10^ — 10^^ M©) at high redshift may have 
very low star formation efficiencies, as a consequence of lower 
metallicities (e.g., Gnedin & Kravtsov 2010, Krumholz & 
Dekel 2011), H2 self-regulation effects (e.g., Kuhlen et al. 
2012), or reduced cooling accretion rates due to local sources 
(e.g., Cantalupo 2010). If their star-formation efficiencies are 
suppressed, these "dark" galaxies can maintain a reservoir 
of gas that will be available for later star-formation in more 
massive systems. This phase of suspended activity may be 
required in order to simultaneously account for several prop- 
erties of the current galaxy population, like the observed 
slopes of the SFR-mass and the TuUy-Fisher relation, as 
has been suggested by Bouche et al. (2010). 

The potential importance of this early phase of galaxy 
formation raises the obvious question of how to discover 
and study such gas clouds that do not have significant star- 
formation. One possibility is try to detect them in HI absorp- 
tion against background quasars (see, e.g.. Ranch 1998 and 
Wolfe, Gawiser & Prochaska 2005 for reviews) . For instance, 
recent studies by Fumagalli, O'Meara & Prochaska (2011) 
and Cooke et al. (2011) have revealed high column-density 
systems that have properties which may be compatible with 
these proto-galactic clouds, given their metallicity consis- 
tent with primordial or Population III stars. Unfortunately, 
these one-dimensional studies cannot directly discriminate 
between truly isolated clouds or components of larger galac- 
tic reservoirs, either around or within galaxies. Detecting 
these systems in emission rather than in absorption would 
provide crucial information about their sizes, morphologies 
and masses, potentially leading to a direct constraint on 
their star formation efficiencies. 

Direct imaging of HI in the 21cm line is unfortunately 
restricted to the local Universe (e.g. Giovanelli et al. 2005), 
given the intrinsic faintness of this atomic transition. Obser- 
vational attempts to detect "dark" galaxies with this tech- 
nique in the local Universe have not so far produced posi- 
tive results (e.g. Gavazzi et al. 2008). A significant fraction 
of local dwarf galaxies do show, however, a much lower star 
formation efficiency than is found in more massive spiral 
galaxies (Geha et al. 2006, Schiminovich et al. 2010). Ex- 
tending these direct imaging studies to high-redshift, where 
the cosmic star formation rate peaks, may require however 
an intrinsically brighter tracer of cosmic hydrogen: i.e. the 
HI Lya line. 

It has been predicted for many years (e.g., Hogan & 
Weymann 1987) that the neutral gas that is responsible for 
quasar absorption lines at high redshift should produce po- 
tentially detectable fluorescent Lya emission at 1216A from 
the incident ionizing radiation field at A < 912A. For re- 
gions that are optically thick to ionizing radiation and are 
isotropically illuminated from outside, a thin layer around 



the outer surface should act as a simple fluorescent mirror, 
converting up to 60% of the incident ionizing radiation into 
Lya photons (Gould & Weinberg 1996). Velocity field and 
geometrical effects can reduce the signal by half (Cantalupo 
et al. 2005). Unfortunately, the cosmic UV background is 
rather weak - whether predicted theoretically (e.g., Haardt 
& Madau 2012) or estimated observationally using the prox- 
imity effect in quasar absorption spectra (e.g., Bajtlik, Dun- 
can & Ostriker 1988, Ranch et al. 1997, Scott et al., 2002, 
Bolton et al. 2005, DaU'Aglio et al. 2008, D'Odorico et al. 
2008). This makes the detection of fluorescent emission from 
optically thick clouds in the field, i.e. illuminated by the 
general cosmic UV background, extremely challenging, and 
perhaps impossible, even on 8-10 m class telescopes (see, 
e.g.. Ranch et al. 2008). However, local enhancements in the 
ionizing background can be used to substantially increase 
the signal (e.g., Haiman & Rees 2001). For instance, gas in 
the vicinity of a bright quasar may be exposed to a stronger 
UV flux compared with to an "average" cloud, and are then 
expected to be correspondingly brighter in fluorescent Lya 
(e.g., Cantalupo et al 2005). 

The detection of the line-of-sight proximity effect in 
quasar spectra provides evidence - in absorption - that bright 
quasars are indeed "illuminating" and ionizing gas clouds 
within several comoving Mpc scales, at least along our line- 
of-sight. Searches for a similar effect along the transverse di- 
rection using projected pair of quasars, the so called "trans- 
verse proximity-effect", have however produced mixed re- 
sults so far (e.g., Jakobsen et al. 2003, Hennawi & Prochaska 

2007, Worseck et al. 2007, Goncalves, Steidel & Pettini 

2008, Kirkman & Tytler 2008), indicating possible emission 
anisotropics or short ages for the quasar bright phase (< 1 
Myr). Detection of quasar fluorescence - the "counterpart" 
in emission of the proximity effect - may thus provide direct 
constraints on angular emission properties and lifetimes of 
these sources. 

1.1 Previous observational attempts to detect 
fluorescent emission near quasars 

In the last years, we have seen an increase of our knowl- 
edge of the quasar fluorescence process with detailed nu- 
merical studies (e.g., Cantalupo et al. 2005, KoUmeier et al. 
2010). However, until very recently, dedicated observational 
surveys have been limited to a few, generally unsuccessful, 
searches. The first fluorescent survey around a quasar, based 
on the simple expectations from the slab model of Gould 
& Weinberg (1996) was performed by Francis & Bland- 
Hawthorn (2004) using a tunable filter imager on the 4m 
Anglo- Australian Telescope. This produced a null result. 
However, as demonstrated by Cantalupo et al. 2005 (see 
also the discussion in Cantalupo et al. 2007), this null result 
is in fact quite consistent with theoretical expectations from 
quasar illumination considering: i) a more realistic density 
distribution than simple slab geometry for the gas, ii) the 
proximity effect of the quasar itself, which reduces the sizes 
and number density of optically thick cloud^. A few years 
later, Adelberger et al. (2006) discovered, serendipitously, 

^ These effects explained the null detection of fluorescent sources. 
However, it is still unclear why no line emitting galaxies were 
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a double peaked Lya emitter associated with a Damped 
Lya absorbing system in the proximity of a bright quasar. 
The Lya profile was consistent with expectations from flu- 
orescence and indicated that the gas was certainly optically 
thick. However, this objects was also associated with a star 
forming galaxy that was bright enough in the continuum 
to entirely power - internally - the observed Lya emission, 
given its rest-frame Equivalent Width (EWo) of about 75 A, 
well within the range that can be produced by conventional 
stellar populations (e.g., Schaerer 2002). A systematic sur- 
vey for other similar systems in quasar pairs by Hennawi 
et al. (2006) has until now yielded only one Lya emitting 
source (Hennawi et al. 2009). However, these authors con- 
cluded that this object is likely intrinsic to the quasar host 
and similar to the Lya fuzz detected around many active 
galactic nuclei (e.g., Heckman et al. 1991; see also Francis 
& McDonnell 2006). 

Motivated by our own more sophisticated theoretical 
modeling, in Cantalupo et al. (2007) we carried out a blind 
spectroscopic survey around a bright quasar at z = 3.1 us- 
ing the VLT/F0RS2 instrument in a "multi-slit plus filter" 
mode (Crampton & Lilly 1999), using two OIH narrow band 
filters. This method allowed us to perform a deep spectro- 
scopic survey down to a flux limit of Fiya ~ 2 x 10~^* erg s^^ 
cm~^ over a significant but non-contiguous volume around 
the quasar (~ 1700 comoving Mpc^). In good agreement 
with theoretical expectations, we found 13 line emitters con- 
sistent with Lya, a third of which showed signatures of the 
double-peaked profiles expected from fluorescence. Only 2 of 
these 13 objects had a significant {2a) detection of under- 
lying stellar continuum emission in a 2 hour deep V-band 
image of the same field. Unfortunately, the depth of the V- 
band image was not sufficient to rule out the presence of 
internal "illumination" by star formation as the source of 
the ionizing radiation. However, other indications, including 
the constraints on the Lya surface brightness as a function 
of distance from the quasar, suggested that about half of the 
sample could have been indeed fluorescently illuminated by 
the quasar. 

A blind spectroscopic survey has the advantage of re- 
ducing the sky-background and giving immediate confirma- 
tion of the line-emitting nature of the source. However, there 
are also several limitations of this approach, including: i) a 
reduction in the sampled volume, by a factor ~ 14 for our 
"multi-slit plus filter" configuration, compared to narrow- 
band imaging on the same instrument, ii) an inability to de- 
tect extended objects, e.g. filamentary structure, with sizes 
larger than the slit width (2" in our case), iii) the absence of 
unambiguous information on the sizes and positions of ob- 
jects (since one spatial dimension is effectively lost), which 
are crucial to better constrain the EWo by comparison with 
a broad-band image. 

In this paper, we present a new survey that is based on 
deep narrow-band (NB) imaging on VLT-FORS that over- 
comes these limitations and provides the first large statisti- 
cal sample of cosmic gas clouds that have been fluorescently 
illuminated by a quasar. 

The layout of the paper is as follows. In Section [2l we 



detected in the same field (when at least 6 Lya emitters were 
expected). 



present the survey design, observations, data reduction and 
Lya candidate selection. In Section [3l we show the observa- 
tional results. In Sectional we demonstrate that many of the 
objects detected have likely been fluorescently illuminated 
by the quasar. In particular, we compare the EW distri- 
bution and Lya luminosity function with those from field 
surveys, and examine the relation between object luminos- 
ity and the projected distance from the quasar. In Section 
[S] we discuss the nature and physical properties of this new 
population of Lya sources and the implications of our identi- 
fication as fluorescent emission. In Section [6l we discuss the 
uncertainties and limitations of our survey. We summarize 
and conclude the paper in Section [T] Through the paper, we 
use a "standard" ACDM cosmology with Ho = 71 km s~^ 
Mpc"\ flm = 0.27 and Qa = 0.73 from WMAP seven-year 
release (Komatsu et al. 2011). 

2 THE DATA 

2.1 Survey Design 

Our survey design has been motivated by two simple but 
crucial requirements: i) to maximize the fluorescently illu- 
minated volume, by selecting one of the brightest quasars in 
the sky, ii) to minimize the sky-background by having the 
narrowest bandwidth fllter that contains such a volume. In 
order to meet these two requirements, it is vital that the 
systemic redshift of the quasar is known with good accu- 
racy, e.g. from low-ionization lines or narrow forbidden lines 
such as [OIH]. Unfortunately, not many quasars have such 
a measure - and our previously observed quasar at z = 3.1 
is not one of these - especially if we also require them to be 
within an existing NB fllter. Although we had already sev- 
eral detections around the z — 3.1 quasar, a large fraction 
of these would have been too faint to be re-detected in NB 
imaging. For the above reasons, we designed a new custom 
NB fllter to take advantage of the new blue sensitive E2V 
detector on FORS and enabling us to select an ultralumi- 
nous quasar at z = 2.4. Due to the lower surface brightness 
redshift dimming - a factor of about 2 - and the increased 
throughput, an imaging survey based on a 20 hour NB im- 
age a.t z — 2.4 should reach a similar depth to our previous 
spectroscopic survey aX z = 3.1. As we will show, the much 
larger volume probed by NB imaging results in a substan- 
tially larger sample of Lya sources - about a hundred - allow- 
ing us to perform crucial statistical analyses to confirm the 
fluorescent nature of these sources, e.g. the distribution of 
Lya EW compared with similar surveys that did not target 
bright quasars. Moreover, we are now able to constrain the 
sizes and masses of these systems and to detect extended, 
circumgalactic emission. 

2.2 Observations and data reduction 

Observations were taken during four visitor-mode nights at 
the VLT 8.2m telescope Antu (UTl) on 2010 September 9- 
12 using the FORS instrument in imaging mode with the 
blue sensitive E2V CCD. A custom-built interference fllter 
(FILT_414_4) with central wavelength Anb = 414.5nm and 
full-width-half-maximum FWHM= 4nm has been used to 
image the field of HE0109-3518 {zsys = 2.4057 ± 0.0003, 
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bj = 16.7), sensitive to HI Lya over a narrow (A2: = 0.033) 
redshift range around the quasar. The quasar systemic red- 
shift is precisely constrained by the wavelength of the de- 
tected [OIII] A5007 emission line (Shemmer et al. 2004, 
Marziani et al. 2009). The quasar is radio-quiet. The field 
was also observed in the Bhigh and Vhigh filters. A total of 
20 hours integration through the narrow-band filter, 6 hours 
through the V-band (of which 4.5 hours were used) and 0.9 
hours through the B-band filter were obtained. The observa- 
tion log is reported in Table 1. Observations were split into 
separate exposures of 1200 s in the narrow-band and 120- 
300 s in the broad bands. Individual exposures were shifted 
and rotated with respect to each other using a semi-random 
shifting pattern of radius 15" and four different position an- 
gles (0, 90, 180 and 270 degrees) in order to facilitate the 
removal of cosmic rays, CCD cosmetic, ghosts and residual 
flat-field errors. The first part of the night beginning 2010 
September 12 suffered from poor seeing conditions, affect- 
ing 18 V-band exposures (for a total of 5400s). These have 
been discarded before image combination. Therefore, the fi- 
nal V image corresponds to a combined integration time of 
4.5 hours. 

The images have been reduced using standard routines 
within the reduction software IRAF, including bias subtrac- 
tion, fiat fielding, CCD illumination correction and image 
combination. A combination of twilight sky flats and unreg- 
istered science frames has been used to produce flat field 
images and illumination corrections for each band and posi- 
tion angle. Sky subtraction was performed, image by image, 
fitting a polynomial surface function of order two in order to 
compensate for small residuals and possible sky variation. 
Before image combination, each image was registered and 
corrected for distortion and rotation with the IRAF tasks 
"geomap" and "geotran" using a series of unsaturated stars 
and a polynomial surface fit of order 3. Finally, for each 
band, the corrected frames were combined with an averaged 
sigma-clipping algorithm. 

The combined science images have been registered on 
the ICRF frame of reference using the USNO-Bl catalogue. 
The uncertainty of 0".2 in the USNO-Bl catalogue domi- 
nates the astrometric error, except at the edge of the field 
of view, where residuals from distortion correction increase 
the error to about 0".5. The photometric calibration was 
performed using several spectrophotometric standard stars 
from the catalogue of Oke (1990). For the calibration of the 
NB image we have chosen three standard stars (LDS749B, 
LTT9491, NGC7293) that did not present any spectral fea- 
ture in the relevant wavelength range. The derived zero 
points were consistent with each other within few percent. 
Magnitudes quoted here are AB magnitudes corrected for 
galactic extinction (0.09 magnitudes for NB and B band 
and 0.07 magnitudes for the V band; Schlegel et al. 1998). 
Given the large number of exposures of the same field at dif- 
ferent air-masses we were able to determine the atmospheric 
extinction coefficients for each night using relatively bright, 
unsaturated stars. 



Table 1: HE0109-3518 field observation log 



Date 


Piltcr 


Seeing*^ 


4- 

''cxp \ 


2010 Sep. 9 


Vhigh 


0".6 


5400 




FILT_414_4 


0".7 


19500 




Bhigh 


0".6 


720 


2010 Sep.lO 


Vhigh 


0".7 


5400 




FILT_414_4 


0".7 


21600 




Bhigh 


0" .7 


840 


zUiU oep.ii 


Vhigh 


n" "7 
U . i 


ODUU 




T?TT /1 1 /I A 




ii I oU 




Bhigh 


0".6 


480 


2010 Sep.l2 


Vhigh 


1".370".8 


7200 




FILT_414_4 


0".8 


19200 




Bhigh 


0".7 


1200 


TOTAL 


FILT_414_4 


0".7 


72080 




Vhigh 


0".7'* 


21600/16200'' 




Bhigh 


0".7 


3240 



" averaged value. 

* total exposure time per niglit and filter, 
average value for the first exposures up to 5400s. 
discarding bad seeing (> 1") exposures. 



2.3 Lya candidate selection 

The selection and photometry of the Lya candidates have 
been performed using an iterative method based on the pro- 
gram SExtractor (Berlin & Arnouts 1996; v2.5) and IDL 
procedures. The narrow-band (NB) image alone has been 
used for candidate detection. Aperture photometry was sub- 
sequently performed on both NB and broad-band images 
using SExtractor dual-mode. This strategy was preferred 
above a combination of narrow and broad-band as detection 
image (as sometimes done in the literature for Lya emitter 
searches, e.g., Fynbo et al. 2003) since this would bias our 
results towards objects with continuum-detection, while we 
are especially interested in detecting object without signifi- 
cant associated continuum. 

In order to minimize spurious detections we varied the 
main detection parameters in SExtractor creating a large 
number of candidate catalogues. In particular, we used a 
minimum area for detection (DETECT_MINAREA) rang- 
ing from 4 to 8 pixels and a detection threshold (DE- 
TECT.THRESH) spanning from 1.3 to 2 times the local 
noise level. These parameter sets were applied to both un- 
smoothed images and to images smoothed with two different 
gaussian kernels. Each SExtractor pass produced, on aver- 
age, about 3000 detections. The same procedure was applied 
to detecting candidates in the broad-band images. These 
were associated by SExtractor, whenever possible, to the ob- 
jects found in the narrow-band using the dual-image mode. 
Initial photometry was performed with isophotal apertures 
(FLUXJSO) and SExtractor FLUX_AUTO for each band. 

As a first step in the selection of Lya candidates, we 
derived an initial guess on the candidate observed EWs from 
the magnitudes obtained by SExtractor: 



EW = AAnb 



/3a + 2 



X (10 



-0.4(NBab-Vab) 



1) 



(1) 



Where AAnb and Anb are, respectively, the NB filter 
FWHM and central wavelength. Ay is the Vhigh filter effec- 
tive central wavelength, fix is the slope of the UV continuum 
(in units of erg s~^ cm~^ A) between Av and Anb, NBab and 
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Vab are, respectively, the NB and V magnitude obtained by 
Sextractor. We used the V-band to estimate the continuum 
since it is not affected by the presence of the Lya emission 
line. With this strategy we are able to better select objects 
without significant continuum emission. At the same time, 
this simplifies the calculation of the EW. When both are 
present, we used the V-band and B-band isophotal magni- 
tude to derive /3a, otherwise we assigned to the source a flat 
continuum slope (in frequency), i.e. I3\ — —2. This is the 
typical maximal value for a starburst galaxy (e.g., Meurer 

et al. 1999). Lya surveys at 2.1 < z < 2.3 find fix 1.9 

for the most luminous sources and higher values (/3a '^-^ —1.5 
but with a large scatter) for the faintest (e.g., Guaita et al. 
2011, Nilsson et al. 2009). Detected objects with /3a < -3 
are extremely rare. Note that a value of /3a = —3 would 
only imply a decrease of about 25% of our estimated EW. 
For the objects without broad-band detection we assigned an 
initial continuum flux corresponding to the average 2a level 
within a 1" photometric aperture across the image (larger 
than the typical seeing of 0".7). To convert the NB mag- 
nitude to a line fiux, we assumed that the emission line is 
fully contained by the NB filter and centered at the filter 
transmission pealjj. 

Since the filter transmission curve has a roughly gaus- 
sian shape, the derived flux is thus a lower limit to the actual 
emission. This implies that the quoted EW should also be 
considered as a lower limit, strengthening the results dis- 
cussed in section 3. We did not attempt to remove the line 
contribution from the B-band magnitude used to estimate 
/3a . This is also a conservative assumption for the estimation 
of the EW. 

For each catalogue we then performed a first cut se- 
lecting objects with observed EW> 68A, corresponding to 
a rest-frame EWo = EW/(1 + z) > 20A for Lya at z = 2.4. 
We chose this EWo cut in order to be comparable with the 
majority of Lya surveys in the literature. This produced 
Lya candidate catalogues of which the largest contained 224 
objects. Then we cleaned through visual inspection the spu- 
rious detections due to the vicinity of bright stars, exces- 
sive fragmenting of objects (e.g., galaxies showing multiple 
spatial components in the narrow-band and a single spatial 
component in the broad-band) and noise at the edges of 
the image. This produced a set of catalogues containing a 
number of candidates varying from 89 to 100. 

For the objects in the cleaned catalogues without Sex- 
tractor broad-band detections we performed again photom- 
etry using our own specifically designed IDL procedures in 
order to derive an upper limit on the broad-band flux and 
therefore a lower limit on the EW. The sources without Sex- 
tractor broad-band detection are typically spatially compact 
or unresolved in the NB image. Therefore we used circular 
apertures with a flxed size of 1.5 arcsec (twice the seeing 
disk, 0.7", that is coincidentally the same for each band) for 
both the NB and broad-bands. Using exactly the same aper- 
ture for NB and broad-band is again a conservative assump- 
tion for the calculation of the EW, since the gas extended 



^ This is a similar approach with respect to the one used by 
other deep surveys not targeting a bright quasar (e.g., Grove et 
al. 2009, Hayes et al. 2010) that we will use in section 4 for a 
direct comparison of the Lya luminosity functions. 
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Figure 1. Rest frame equivalent width (EWq) of the detected ob- 
jects versus their narrow-band magnitudes. Blue solid (or open) 
circles represent candidates with V-band (or both V and B-band) 
detections above the 3cr level for V-band (and 5cr for the B-band). 
When the B-band is detected (open circles) the EWo is derived 
from a continuum conservatively estimated from both V and B- 
bands without correcting for the line contribution to the B-band. 
For objects without B detection, we estimate the continuum from 
the V-band, assuming a continuum slope of /3a = —2. The red 
solid circles represent lower limits on the EWo of the objects 
without continuum detection (i.e., V< 3a"). These correspond to 
the (local) Icr lower limits for objects with signal to noise ratio 
SNR< in V-band (the deepest available), otherwise the mea- 
sured flux plus the local Icr noise has been used (e.g., Feldman 
& Cousins 1998; see text for details). The dashed line represents 
the theoretical limit on the EWo that is produced by "normal" 
stellar populations (EWo = 240A ; e.g.. Chariot & Fall 1993, 
Schaerer 2002) while the dotted line is the lowest EWq for candi- 
date selection (EWo = 20A , consistent with most published Lya 
surveys). The continuum-undetected objects with EWo > 240A 
(red dots above the dashed line) are the best candidates for proto- 
galactic, or "dark" , gas-rich clouds that have been fiuorescently 
illuminated by the quasar. 



distribution and scattering effects would likely result in a 
Lya emission that is more extended than the stellar contin- 
uum (e.g., Cantalupo et al. 2005; see also Steidel et al. 2011). 
The apertures were centered, for all bands, at the location 
corresponding to the peak of the NB emission (derived from 
SExtractor) . We used relative astrometry to refine the posi- 
tion of these apertures in the broad-bands. If the flux mea- 
sured in the broad-band apertures was above 3 times the 
local noise level (a) we considered the source continuum- 
detected and we used the measured flux to compute the 
EW. Otherwise, we consider the object undetected in the 
continuum and we estimate an upper limit on the contin- 
uum fiux using a procedure similar to Feldman & Cousins 
(1998). For apertures where the measured flux is positive 
this corresponds to assigning an upper limit equal to the 
sum of the measured value and the local la. In case the flux 
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is negative, although Feldman & Cousins (1998) provide a 
table for assigning upper limits (with values < la) in this 
situation, we decided to use as a conservative value the local 
la as a minimum threshold. The continuum flux upper lim- 
its translates then into EW lower limits. With this method, 
we use the statistical information contained in the measured 
flux and we have a consistent measurement with the la er- 
rorbars of the detected objects, i.e. we make sure that the 
actual EW is above the reported value within at least la 
confidence level. We note that current Lya surveys in lit- 
erature use different methods to estimate continuum flux 
upper limits for undetected sources. In particular, a com- 
mon choice is to use a flx value (e.g., a globally estimated 
la or 2a value for a given aperture) independently of the 
measured flux. This complicates the comparison of different 
works in the literature. 

Bright foreground objects in proximity to the aper- 
ture were masked, but we cannot rule out the possibility 
that fainter foreground objects are contaminating the aper- 
tures in the broad-band. Note that this would imply that 
the derived EW of the objects in question are underes- 
timated and thereby strengthen the conclusions presented 
in the paper. For the objects with SExtractor broad-band 
detection we used FLUX_AUTO to calculate total magni- 
tudes and isophotal apertures (FLUXJSO) to compute the 
colors and EW, as commonly performed in previous Lya 
searches (e.g., Hayes et al. 2010). Isophotal apertures and 
FLUX_AUTO measurement did not produce magnitudes 
varying more than few percent from each other, except 
in the few cases where sources were located close to fore- 
ground objects. Only objects with a SNR>5 in the narrow- 
band images were retained in the catalogues. This corre- 
sponds to about NBab(5(7) ~ 26.8 on average across the 
image, while the typical (median) la magnitude limits for 
V-band and B-band were, respectively, Vab(1o") ~ 29.7 and 
BAB(lcr) ~ 29.1 for a 1.5" diameter circular aperture. 

In order to assess the level of spurious detections due 
to noise for each set of parameters and filters we produced a 
negative catalogue applying SExtractor to the narrow-band 
image multiplied by —1. The number of detected objects 
in the negative catalogues varied from zero to 30. In the 
latter case, however, these spurious detections were mostly 
confined at the image edges, while only between three to 
five spurious objects were present at the center. Among the 
eleven candidates that were not present in every positive 
catalogue, nine did not have any broad-band detection. For 
these nine candidates we have rechecked the narrow-band 
image splitting the total exposure time in four parts (each 
corresponding to a given rotation angle). Six objects were 
detectable (at SNR>3) in at least three out of four images. 
The remaining three objects, in accordance with the number 
found in the negative catalogues, were considered spurious 
and thus removed from their catalogues. The final catalogue 
contained thus 97 sources. In section IGTTI we discuss in detail 
the possible fraction of foreground contaminations in our 
sample. 

Finally, we estimated the completeness of our NB imag- 
ing with the standard procedure of distributing and recover- 
ing a large number (~ 1000) of artificial sources with varying 
sizes. In particular, we used the distribution of apertures 
derived by SExactror for our detected candidates obtain- 
ing a completeness factor fc = (25%, 50%, 75%, 90%) for 



NB = (26.95,26.55,26.30,26.12) corresponding to Lya lu- 
minosities, at z = 2.4, LLya = (0.20, 0.30, 0.38, 0.43) in units 
of lO'*^ erg s^^ . The factor fc will be used to correct the Lya 
luminosity function for incompleteness in section |4?2] 



3 RESULTS 

The final, merged catalogue consists of 97 narrow-band ex- 
cess objects. The position and full properties of all detected 
objects are reported in tabular form in Tables O [3] and 
Appendix A. The EW (estimated as in section 12.31 ) and 
NB magnitudes of these objects are presented in Figure [1] 
Among the 97 NB-detected objects, 31 do not have any de- 
tectable continuum (above 3a) in the very deep broad-band 
images (red symbols with arrows in Figure [1} and there- 
fore the reported EW has to be considered a lower limit as 
discussed in section [2^ Given the sensitivity of our broad- 
band imaging we are able to put strong constraints on the 
EW limits, at least for the objects with NB magnitudes 
NBab < 26. Among these, 12 Lya candidates have a lower 
limit EWo > 240A, i.e. higher than the typical maximum 
value produced by normal stellar populations (e.g.. Char- 
lot & Fall 1993, Schaerer 2002; see also Malhotra & Rhoads 
2002 and Raiter, Schaerer & Fosbury 2010). Six of the 66 
objects that have at least one broad-band detection (blue 
circles in Figure [l]) also have EWo > 240 A. As we will show 
in this section, these contain also the objects with the most 
evident signs of extended, filamentary emission. 

In Figure (2] we present the combined narrow-band im- 
age (20 hours combined integration time) with overlaid posi- 
tions of the 97 Lya candidates. The quasar is near the center 
of the image. Open squares (circles) represent objects with a 
measured or lower limit EWo > 240A (EWo < 240A) while 
red (or blue) colors indicate if the object was undetected (or 
detected) in any of the broad-band images. 

From the image, we have visually identified an addi- 
tional narrow-band excess source, i.e., a blob-like extended 
emission (magenta ellipse in Figure [2] approximatively at 
RA = 1 : 11 : 57 and S — —35 : 01 : 50). Due to the presence 
of relatively bright components within this extended image, 
which appear also in the broad-band (possibly, foreground 
objects), this extended source was split into several parts 
by SExtractor and it did not make it to the final catalogue 
using our selection technique. We will refer to this object as 
LABI from now on. 

In Figure [3] we show postage stamp images of 12 can- 
didates without continuum detection {V < 3a) and lower 
limit EWo > 240 A. The objects are ordered by decreasing 
EW limit from the top to the bottom panels. Their posi- 
tions, photometric properties and EW limits are reported in 
Table 13 Among the brightest candidates, #1484 and #902 
have a very high lower limit on their rest-frame Lya equiva- 
lent widths, namely EWo >> 400A. Interestingly, #1484 is 
also the source closest to the quasar (in projection), located 
just about 15 arcsec away (the quasar is just above the top 
right corner in the first panel of Figure 3). Note that the 
vast majority of the objects without broad-band detection 
appear spatially compact. 

In Figure 21 we show a stacked NB (left-hand panel) 
and V-band (right-hand panel) image of the 12 proto- 
galactic candidates, i.e. the continuum-undetected sources 
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Right ascension 



Figure 2. The combined narrow-band image (20 hours total integration time) of the field surrounding HE0109-3518 (near the center of 
the image). Squares (or circles) indicate the position of detected emission line objects with EWo > 240A (or EWo < 240A), assuming 
the line is hya. Red (or blue) colors represent objects undetected (or detected) in the broad-band images. The position of the extended, 
blob-like emission source (LABI) is indicated by a magenta ellipse (approximatively at ra = 1 : 11 : 57 and dec = —35 : 01 : 50). Objects 
with undetected continuum and with EWo > 240A (i.e. the red squares) are the best candidates for protogalactic clouds, or "dark" 
galaxies, fluorescently illuminated by the quasar (see text for details). 



with EWo > 240 A. Before stacking, bright sources have been 
masked everywhere except in the central 3" -wide region (in- 
dicated by the black circle) and the local (30" radius), sigma- 
clipped mean has been subtracted from each image. We then 
performed aperture photometry within a 2" diameter aper- 
ture on both the NB and V-band stack, centered on the po- 
sition derived from the NB. The measured flux in V-band is 
below the la level for a 2" aperture, i.e. Vistack ~ 30.3, imply- 
ing a combined constraint of EWq, stack > 800A (la). This 
measurement strongly suggests that internal stellar sources, 



if present at all, do not contribute signiflcantly to the Lya 
emissiorlf]. The stacked Lya emission proflle remains com- 



As a reference, we have estimated with CLOUDY (vlO.OO; Fer- 
land et al. 1998) that an EWq > 800A is still one order of magni- 
tude lower than the value expected from pure nebular continuum 
if measured at rest-frame I6OOA and using = —2 for a fully 
ionized cloud by a typical quasar spectrum {a = —1.7; e.g. Teller 
et al. 2002). 
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Figure 3. Images of the narrow-band excess sources without continuum detection {V < 3cr) and EWq > 240A (assuming the line is 
Lycf) in three different bands: narrow-band (left-hand panels), B (central panels), and V (right-hand panels). For an object at z = 2.4, 
the central wavelength of these filters correspond to the rest-frame Lya, 1300A and I6OOA, respectively. Note that the B-band contains 
also the source Lya line, differently from the V-band, although the emission would be typically too faint to be detected in B. Panel 
dimensions are 15 X 15 arc sec^. The images have been smoothed with a gaussian kernel with radius 0".5 (two pixels) for clarity. These 
objects represent the best candidates for fluorescent protogalactic clouds or "dark" galaxies given the absence of detected continuum 
and the high lower limit on the EWo- As discussed in the text, an EWo > 240A strongly suggests that internal star formation is not 
likely the origin of the Lya emission and that we are detecting dense gas fluorescently illuminated by the quasar. 
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Figure 4. Stacked NB image (left-hand panel) and V-band image (right-hand panel) of the 12 Lya sources that did not have individual 
continuum detection and an implied EWo > 240 A. The images have been smoothed with a gaussian kernel with radius 0".5 (two 
pixels) for clarity. The V-band stack does not show any continuum source above la level (V~ 30.3). The combined constraint implies a 
EWo > 800A, effectively ruling out internal star-formation as the origin of the Lya emission in these objects. 



pact and consistent with a point source with FWHM size of 
about 1" (about 8 physical kpc at 2 = 2.4). 

Postage stamp images and summary properties of the 
Lya candidates with broad-band detections (> 3a) and 
measured EWo > 240A are presented in Figure [S] and Ta- 
ble [31 respectively. It can be seen that the three brightest 
Lya candidates in our sample, #1159, #2886 and #1412 
are also the objects with the highest measured EW, i.e. 
EWo - 350A. They are detected both in V and B bands 
and show a very red V — B color. In particular, their B 
magnitudes are faint enough that the very bright Lya emis- 
sion could contribute most of the broad-band flux. The pres- 
ence of a stellar continuum is favored by the clearly detected 
flux in V-band, but either the stellar population is old or it 
is embedded in a large amount of dust that is not affect- 
ing the Lya emission. This would be necessary to reconcile 
the red color with the high Lya EW. As we will discuss 
in section [531 a possible explanation is that quasar fluores- 
cence is illuminating a large, dust-poor gas reservoir - the 
Circum Galactic Medium (CGM) - surrounding dust-rich, 
star-forming galaxies. 

At least three objects, #1159, #1412 and #2411, show 
evidence of extended, possibly filamentary, emission. In par- 
ticular, the image of #1159 shows a hint of three structures 
extending out of the bright nucleus. The longest of these has 
a (projected) linear size of about 4" , corresponding to about 
30 physical kpc a,t z — 2.4. This scale is similar to the virial 
radius of a ~ 10^^ dark- matter halo at this redshift. 
Object #2411 also shows a filamentary "bridge" on similar 
(projected) scales apparently connecting one region of ex- 
tended emission without a compact nucleus (and with only 
a marginal detection in V-band) to another bright narrow- 
band excess object with high EW (#2433 ). Revealing the 
nature of these two apparently connected Lya companions is 
complicated by the fact that the broad-band images show a 
continuum source between them. Without knowing the red- 
shift and the spectral energy distribution of this continuum 
object, it is difficult to estimate its contribution to the NB 
observed flux. However, the red B — V color and the dif- 
ferent morphology argue that such compact object cannot 
be the only source of narrow-band emission in the "bridge" . 
If #2411 and #2433 are physically connected, this emitting 
"bridge" may have a similar origin to the narrow gas filament 



discovered by Ranch et al. (2011) in a deep blind spectro- 
scopic survey at 2 = 3.3 and corresponding to a continuum 
source showing tidal tails in HST imaging. 

In Figure [6] we show the narrow and broad-band 
postage stamp images of the blob-like extended emitter 
LABI. Morphologically, this object resembles other elon- 
gated Lya blobs previously discovered (e.g., Steidel et al. 
2000, Matsuda et al. 2004, Nilsson et al. 2006, Yang et al. 
2009, Prescott et al. 2012). It shows a relatively uniform 
surface brightness except for one peak, which is off centered 
with respect to the extended emission. At least four contin- 
uum sources are detected in the broad-band image(s) within 
the projected area of LABI. One of these sources is located 
in the proximity of the NB emission peak. The overall shape 
of LABI seems to follow these continuum objects but at 
present we are not able to confirm if they lie at the same 
redshift as LABI. Note that the bright NB excess object in 
the top-left corner in Figure [6l (about 10" away from LABI) 
is #2620, one of the objects without continuum-detection 
with high EW limit (EWo > 276A ). 



4 PHYSICAL ORIGIN OF THE EMISSION 

In this section, we compare our observational results to those 
obtained in similar deep searches for Lya emission in re- 
gions that do not target high-redshift quasars (we refer to 
these surveys as "blank-field" studies) and with the results 
of radiative-transfer simulations. In particular, we will ex- 
amine the EW distributions, Luminosity Functions and the 
relation between (projected) distances form the quasar and 
candidate luminosities. The combination of these statistics 
provide strong evidence for the fluorescent origin of the line 
emission of many of the objects in our sample. 

4.1 EW distribution 

The Lya EW is one of the most important diagnostics to 
distinguish between internal star-formation and other mech- 
anisms that can produce Lya emission (e.g., Schaerer 2002), 
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Table 2: Position and properties of Lyo candidates without continuum detection {V and B < 3cr) 



Id 




RA 


DEC 


NB" 


V?. 


SNR(V)'= 




pd 

Lya 


Ltya 






(J2000) 


(J2000) 


(AB) 


(AB) 




(AB) 


10^-"^^ erg s^^ cm~2 


10*2 erg s-i 


A 


1484 


1: 


:11:44.83 


-35:03:12.1 


oc- qi +0.09 


>29.71 


-1.1 


>28.16 


1.95± 0.16 


n ni 1 n nT 

0.91± 0.07 


>483 


902 


1: 


:11:41.80 


-35:04:43.5 


„_ p„+0.15 
Z0.0(_o j4 


>29.75 


-0.8 


>29.45 


1 CT w 1 n on 

1.54± 0.20 


n "71 1 n nn 

0.71± 0.09 


>443 


1737 


1 


:11:34.45 


-34:59:41.9 


OK nc+0.08 
^0-U»_o.o8 


>29.22 


0.5 


>29.38 


o o 1 n 1 T 


1 1 o 1 n no 

1.12± 0.08 


>325 


390 


]^ 


:11:43.50 


o CT r\cr AO t 

-35:05:48.7 


OK KV+0.14 


>29.75 


-1.5 


>29.41 


1 Art \ n 1 T 

1.4U± 0.17 


n n tr i n no 

0.65± 0.08 


>315 


2944 


]^ 


■11'45 28 


-35:02:14.1 


25.42+«-iJ 


>29.69 


-0.3 


>29.38 


1.77± 0.17 


0.82± 0.08 


>303 


1321 


1' 


11:51 76 


-35:03:42.0 


25.87t«-I 


>29.78 


-0.3 


>29.15 


1.17± 0.17 


0.54± 0.08 


>297 


555 


]^ 


■11'42 91 


-35:05:24.5 


25.74l«:i? 


>29.80 


-0.6 


>29.47 


1.32± 0.14 


0.61± 0.06 


>282 


2620 


]^ 


.1 1 .^7 m 

. ll.O t .U-L 


-35:01:35.9 


25.001H? 


>29.12 


0.8 


>28.62 


2.60± 0.16 


1.21± 0.08 


>276 


2131 


]^ 




-35:00:27.9 


9^ Ko+0.11 


>29.51 


0.3 


>29.54 


1.52± 0.15 


0.71± 0.07 


>264 


2274 


]^ 


.1 1 .on 1 7 


-35:00:47.3 


26.2410:^6 


>29.86 


-1.8 


>29.54 


0.83± 0.11 


0.39± 0.05 


>263 


1519 


^. 


. 1 1 . Qft Q« 


-35:03:17.0 


25.72t[;-l5 


>29.62 


-0.9 


>29.31 


1.34± 0.17 


0.62± 0.08 


>254 


1212 




■11'45 14 


-35:03:57.1 


26.001?,;^° 


>29.49 


0.3 


>27.91 


1.04± 0.09 


0.48± 0.04 


>240 


1498 


1: 


:11:47.81 


or' no "1 o o 

-35:03:13.8 


o /'a^H'~0.1y 

26.07Io ig 


>29.71 


-1.6 


>29.43 


n nrr 1 n "1 /I 

0.97± 0.14 


n /I r' 1 n n^7 

0.45± 0.07 


>218 


2842 


1: 


:11:49.17 


or' no OA n 

-35:02:34.0 


on on— |— O.IT 

26.26Io 15 


>29.68 


-1.5 


>29.39 


n o 1 1 n T o 

U.81± 0.12 


n o o 1 n n r* 

0.38± 0.05 


>216 


1829 


1: 


:11:30.15 


-34:59:51.5 


oc no-|-0'15 

25.98Io 14 


>29.80 


-0.1 


>28.50 


1 ncT 1 n 1 H 

1.05± 0.14 


n Ar\ 1 n r\c 

0.49± 0.06 


>202 


2985 


1: 


:11:46.39 


o no o "1 r' 

-35:02:21.5 


OK nn+0.14 

25.99Io 13 


>29.69 


-0.1 


>29.38 


1 n 1 n 1 o 

1.04± 0.13 


n A r\ \ n nn 

0.49± 0.06 


>174 


1565 


1: 


:11:46.79 


or* no o o "i 

-35:03:22.1 


oc 1 n+0.17 

26-19Io.i6 


>29.71 


-3.0 


>29.44 


n o ^7 1 n 1 o 

0.87± 0.13 


n A f\ \ n nn 

0.40± 0.06 


>174 


1926 


1: 


:11:33.06 


o r* nn n t o 

-35:00:01.2 


^0-»'_0.12 


>29.47 


0.4 


>28.38 


1 1^1 n "1 o 

1.17± 0.13 


f\ ^ A \ n nn 

0.54± 0.06 


>163 


619 


1: 


:11:58.87 


or* n r' 1 z'* o 

-35:05:16.2 


on n^— pU.XO 

26.04Io 14 


>29.26 


0.6 


>29.54 


"1 nn 1 n "1 o 

1.00± 0.13 


n /I n 1 n nn 

0.46± 0.06 


>156 


2034 


1: 


:11:43.31 


o nn o t o 

-35:00:21.3 


on "1 o— |— 0.15 

26.13Io 14 


>29.57 


0.2 


>29.52 


n no 1 n T o 

0.92± 0.12 


n .i o 1 n nn 

0.43± 0.06 


>146 


1706 


1: 


:11:56.37 


OA r' n on 

-34:59:30.7 


oc no-|~0.16 

25.92Io 14 


>29.39 


-0.2 


>28.50 


l.llit 0.15 


n r' o 1 n n^7 

0.52± 0.07 


>140 


85 


1: 


:11:38.64 


or* f\r^ OT 

-35:06:27.6 


on ^70—1— 0.22 

26.78Io ig 


>29.59 


-1.0 


>29.04 


n r'n I n nn 

0.50± 0.09 


n o o 1 n n .1 

0.23± 0.04 


>137 


1279 


1: 


:11:40.40 


or* no A T 

-35:03:45.7 


o/i nc+0.09 
24.96Ioo9 


>28.26 


2.7 


>29.39 


o '^n 1 n 0"i 

2.70± 0.21 


1 OKI n "1 n 

1.25± 0.10 


>120 


1439 


1: 


:11:52.60 


-OtJ.UO.UtJ.i 


ZD.U4_Q 13 


>29.78 


-0.3 


>29.06 


1 no-i- 1 1 


n Af\-\- n n"^ 


>115 


2454 


1: 


:11:59.18 


-35:01:11.8 


26.62j;0ig 


>29.52 


-0.2 


>28.62 


0.58± 0.11 


0.27± 0.05 


>111 


70 


1: 


:11:45.10 


-35:06:28.1 


26.57l«;25 


>28.58 


1.0 


>28.93 


0.61± 0.13 


0.28± 0.06 


>43 


1432 


1: 


:11:47.69 


-35:03:04.5 


26.40j:«:2g 


>28.21 


2.9 


>29.43 


0.72± 0.14 


0.33± 0.06 


>42 


1762 


1: 


:11:53.86 


-34:59:39.6 


26.53i;5-2i 


>28.53 


1.6 


>29.27 


0.64± 0.11 


0.30± 0.05 


>41 


2525 


1: 


:11:36.41 


-35:01:21.3 


26.50+«:20 


>28.41 


2.7 


>29.50 


0.65± 0.11 


0.30± 0.05 


>32 


665 


1: 


:11:51.87 


-35:05:09.6 


26.56l«;i? 


>28.37 


2.9 


>29.53 


0.62± 0.10 


0.29± 0.05 


>32 


1989 


1: 


:11:38.01 


-35:00:25.4 


26 47+0-31 


>27.82 


2.8 


>28.40 


0.67± 0.17 


0.31± 0.08 


>28 



" Total magnitude. 

^ Lower limits - measured value plus the local la or local la if SNR< - derived from a circular aperture of 1".5 diameter. 
V-band Signal to Noise Ratio within a circular aperture of 1".5 diameter. 

Derived from total NB magnitudes assuming that emission line is at the filter central wavelength (lower limits). 
Derived from Fi^y^ assuming z= 2.4. 
^ Rest frame EW limits within a 1".5 diameter circular aperture estimated from Vum assuming an UV continuum slope I3x = —2. 



Table 3: Position and properties of Lya candidates with continuum detection {V > 3a) and EWo > 240 A 



Id 


RA 
(J2000) 


DEC 
(J2000) 


NB 
(AB) 


(AB) 


(AB) 


Lya 

10^-''' erg s^-*^ cm^^ 


Lya 

1042 erg s-i 


EWg 

A 


1159 
2886 
1412 
2411 
2914 
2936 


1:11:56.59 
1:11:46.88 
1:11:45.42 
1:11:33.48 
1:11:33.93 
1:11:53.84 


-35:03:46.9 
-35:02:33.8 
-35:03:05.3 
-35:01:09.8 
-35:02:13.0 
-35:02:18.1 


22.660^04^^ 

23.180^04,4 

22 47+0-04 
^^■^'-0.04 

24 25+0.06 

25 14+0.08 
Z0.1t_Q 08 

23 qn+o os 
zo.yu_Q 05 


24 74+0.04 
^^•'^-0.04 

25 99+0.07 
zo.yy_o 06 
25 45+0.06 
zo.'iO_o 06 

27.81^0;?^ 

28.69+0-^* 
27.34to:24 


26.06toi2 

26.72 j:o;j7 

25.98 tVA 

> 28.73 

> 28.52 

> 27.98 


22.44±0.81 
13.90±0.50 
26.73±0.97 
5.19±0.28 
2.29±0.16 
7.16±0.32 


10.42±0.38 
6.46±0.23 
12.42±0.45 
2.41±0.13 
1.06±0.08 
3.33±0.15 


358±35 
356±48 
327±28 
296±51 
291±77 
262±48 


2433^ 


1:11:33.85 


-35:01:12.4 


90 09+0.09 

08 


25.35«:- 


26.3 lO-.l^ 


12.22±0.97 


5.68±0.45 


218±27 



" See Appendix A for the list of the remaining 61 Lya candidates with continuum detection. 
' Total magnitudes (or lower limits for B-band, see Table 2) . 

Assuming that emission line is at the filter central wavelength. 

Derived from Fi^y^ assuming z= 2.4. 
" Rest frame EW for Lya at z= 2.4. A UV continuum slope fix = —2 is used for objects undetected in B. 
f (Projected) "Companion" of #2411. 
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Figure 5. Images of the objects with continuum detection {V > 3ct) but still with EWo > 240A (assuming the line is hyalpha) in three 
different bands. These sources include the extended objects with filamentary emission (#1159, #1412, and #2411 with its "companion" 
#2433). The high value of the EWo suggests that local star-formation cannot be the only source of the detected radiation. Instead, the 
emission is fully compatible with dense circum-galactic gas fluorescently illuminated by the quasar (see text for details). Panel dimensions 
are 15 x 15 arcsec^. The images have been smoothed with a gaussian kernel with radius 0".5 (two pixels) for clarity. 



Figure 6. Postage stamp images of the extended, blob-like emission LABI in narrow-band (loft-hand panel), B-band (central panel) and 
V-band (right-hand panel). Panel dimensions are 20x20 arc sec'^. The images have been smoothed with a gaussian kernel with radius 
0" .5 (two pixels) for clarity. The compact narrow-band excess object on the top-left corner of the NB image is #2620, one of the brightest 
fluorescent protogalactic or "dark" cloud candidates (EWq > 276A ) and undetected in both B and V-band). 



simply because of the limit to the ratio of the number of 
ionizing photons to continuum photons at the Lya wave- 
length from possible stellar populations. In the proximity 
of a quasar: the Lya emission may be boosted by fluores- 
cent emission, leaving untouched the stellar continuum, if 
present. The net effect is an EW "boost" of every object 
directly illuminated by the quasar by an amount that may 
vary with the object gas mass and its distance from the 
quasar. Therefore, if a significant fraction of our objects are 
indeed fluorescently boosted, our EW distribution should be 
skewed towards higher values with respect to surveys away 
from bright quasars. 

Slightly different observational methods and selection 
techniques complicate the comparison to surveys in the lit- 
erature. Despite these difficulties, the observed EW distri- 



butions from several recent 'blank-fields' surveys at 2.1 < 
z < 2.8 are remarkably similar. In particular, they can all 
be fitted by a single exponential function with an e-folding 
scale length of «; ~ 50A . These surveys include: Grove et 
al. 2009 {z = 2.8, Liy^ > 3 x 10''^ erg 8"^ original data 
from Fynbo et al. 200^ ) , Nilsson et al. 2009 {z = 2.3, 
Liyc > 2 X 10*^ erg s"^), and CiarduUo et al. 2012 (z = 2.1, 



Extracted from a field containing a damped Ly« system at 
2 = 2.85 (Q2138-4427); we computed the EW distribution using 
all photometrically selected candidates in this field (Tables A3 
and A4 in Grove et al. 2009). Note that, although this is not 
a strictly blank field since it contains a DLA, the over density 
of Lya candidates is only a factor 1.1 with respect to the other 
surveys. 
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Figure 7. The cumulative rest frame equivalent width (EWo) distribution for all the detected line emitting objects in our sample (left- 
hand panel; L > 2 X 10*^ erg s~^) and considering only the candidates with Lj^y^ > 10*'^ erg (right-hand panel), always assuming 
the line is Lyo. Objects are split according to their detection (blue dotted histogram) or non-detection (red dashed histogram) in any 
of the broad-bands. Note that the EWq of the objects without continuum detections are lower limits. The solid black histogram shows 
the cumulative EW distribution including all objects and using the lower limits for the continuum-undetected candidates. The results 
from recent surveys targeting 'blank-fields', at 2.1 < z < 2.8 are shown (black dot-dashed line, rescaled to match the total number of 
objects) taken from CiarduUo et al. 2012 (z = 2.1, L^y^ > lO^'^ erg s~^; original data from Guaita et al. 2010), Nilsson et al. 2009 
{z = 2.3, L^y^ > 2 X 10*2 erg s"!), and Grove et al. 2009 {z = 2.8, L^y^ > 3 X IQ-*! erg s^^; original data from Fynbo et al. 2003). 
The observed EW distributions from these surveys can all be fitted by a single exponential function with an e-folding scale length of 
w ~ 50A (see text for details). Our sample presents a clear excess of high EW objects with respect to the 'blank-fields' surveys. This is 
true also considering only the continuum-detected objects or the more luminous part of the sample (right hand panel). This clear excess 
of high EWo with respect to the field surveys is a distinctive sign of a fluorescent hya boosting due to the quasar proximity. 



Liya > 10*^ erg s^^; original data from Guaita et al. 2010) 
. As shown by these studies, the EW distribution seems 
to be also largely independent of the luminosity cut except 
at the very bright end, where there is an apparent lack of 
high EW objects. The EW distribution appears instead to 
evolve with redshift: at z ~ 3.1 there are evidences for an 
increase in EW scale length to w ~ 70A, (CiarduUo et al. 
2012; Nilsson et al 2009; see also Ouchi et al. 2008) or larger 
at z > 4 (e.g., Malhotra & Rhoads 2002, Dawson et al. 2004, 
Saito et al. 2006, Shimasaku et al. 2006). 

In Figure [T] we present the cumulative EW distribu- 
tion for all the detected objects in our sample (left-hand 
panel; L > 2 x 10*^ erg s~^) and considering only the can- 
didates with Liya > W^^ erg s~^ (right-hand panel). The 
sample is divided according to whether the continuum was 
detected (blue dotted histogram) or not (red dashed his- 
togram). The EWo of the objects without continuum detec- 
tions are of course therefore lower limits. The solid black 
histogram shows the cumulative EW distribution including 
all objects and using the lower limits for the continuum- 
undetected candidates. 

Our sample contains a clear excess of high EW objects 
with respect to the 'blank-field' surveys. This remains true 
also considering only the continuum-detected objects (that 

^ The only exception is given by the results of Hayes et al. (z = 
2.2, Liyct > 3 X 10'*^) who find a slightly larger e-folding scale 
with value w ~ 76A but they suggest that poor statistics are 
likely responsible for this discrepancy. 



show an e-folding scale of w ~ 85A ) or when we consider 
only the most luminous part of the sample: in this case the 
e-folding scale for the continuum-detected objects (the large 
majority) is of the order of ui ~ 120A, although the statistics 
are poorer. In both cases, a K-S test confirms that our EW 
distribution is different from the one obtained by the "blank- 
field" surveys at a significance greater than 99% . It should 
be noted that we that we have not attempted to correct for 
the non-square bandpass of our filter since we have assumed 
that the Lya emission of each candidate lies at the peak 
transmission wavelength of the filter (since we do not know 
the redshift distribution for all of our objects). We stress 
that this is a conservative assumption for the calculation of 
the EW and likely gives an underestimate of the true Lya 
flux and EW of our candidates. 

In addition to explaining the excess of high EW ob- 
jects, the fluorescent Lya boosting due to the proximity of 
the bright quasar is also able to account for the w increase 
with luminosity that is exhibited by those objects with de- 
tected continua, the opposite trend to that exhibited by the 
"blank-field" surveys. Quasar fluorescence boosts both the 
Lya luminosity and the EW of a given object (independent 
of its "original" position in the Ljj^a-EW diagram ) creating 
a positive correlation between these two quantities. 

Is there any other effect that could explain the differ- 
ent EW distribution obtained in our sample with respect 
to the "blank-field" surveys? Differently from the "blank- 
fields", our sample could be associated with a high density 
peak in the matter distribution due to the presence of the 
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Figure 8. The luminosity function (number density per decade 
in luminosity) of our hya candidates compared with other two 
deep surveys in "blank- fields" (see labels within figure). Solid cir- 
cles (or open squares) show the observed Luminosity Function 
including (or excluding) the correction for completeness. The er- 
ror bars show 84% confidence levels based on Poisson and bino- 
mial statistics (Gehrels f986). Note that, consistently with the 
other surveys, we have not corrected for the non-square bandpass 
shape of the NB filter (see text for details). We detect a clear ex- 
cess of sources at the faint end with respect to "blank-fields" . As 
discussed in section 4.2, this is fully compatible with fluorescent 
boosting due to quasar proximity. 



quasar. We discuss extensively this possibility in the follow- 
ing sections. Venemans et al. 2007 performed a survey for 
Lya emitters around high-redshift radio galaxies detecting 
a significant enhancement with respect to the field (about a 
factor 2-5). However, their EW distributions do not appear 
significantly different from the one obtained by the "blank- 
field" surveys. In particular, the fraction of candidates with 
EWo > 240A (and Liy^ > 10"*^ erg s"^) is less than 5% (B. 
Venemans, private communication), similarly to the "blank- 
field" surveys at 2 < 2; < 3. As a reference, we find instead 
that more than 20% of our sources with Liya > lO**^ erg s~^ 
have EWo > 240A . 



4.2 Luminosity Function 

Similarly to the EW distribution, we expect fiuorescent 
boosting to also impact the Lya Luminosity Function (LF) 
of our sample. In Figure [51 we show the LF obtained from 
all the Lya candidates (open squares) and including a cor- 
rection for completeness (solid circles; see section r2.3|) . We 
compare our LF with the results of Grove et al. 2009 and 
Hayes et al. 2010. These two surveys, the deepest avail- 
able at 2 < 2 < 3, reached a comparable depth and used 
a very similar NB filter on VLT-FORS with respect to 
our study. Despite the small volumes sampled, these two 
studies derived very similar luminosity functions, consistent 
with a Schechter function with log(L*/erg s~^)~ 43.1, slope 
a ~ —1.5 and normalization log((/i*/cMpc~'^) ~ —3.9. 

In our own survey, we detect a clear excess of Lya can- 
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Figure 9. Cumulative luminosity function (CLE) of our Lyo can- 
didates (black solid histogram) compared with other deep surveys 
in the field (see labels within figure). Blue dotted (or red dashed) 
histogram shows the CLF including only objects with (or with- 
out) continuum-detection. 
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Figure 10. Simulated cumulative luminosity function (CLF) ob- 
tained from simulations assuming that the quasar is off (blue 
dotted histogram), i.e. no increase in the UV background due to 
quasar proximity, and that the quasar is turned on (red dashed 
histogram). These simulations include full three-dimensional ra- 
diative transfer of ionizing and Lya photons (Cantalupo & Por- 
ciani 2011; Cantalupo et al. 2005) within high-resolution cos- 
mological simulations of structure formation around a massive 
halo (see text for details). The overall vertical normalization of 
the CLF of Lyo emitters in the QSO-OFF simulation has been 
rescaled to match the observed field CLF of deep surveys (Grove 
et al. 2009, Hayes et al. 2010), while the slope of the CLF is repro- 
duced naturally by our simulation. Using the same normalization 
factor the numerical model is able to reproduce the increase in 
the total number of objects and the steepening of the CLF as 
observed in our data once the quasar is turned on. 
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didates with respect to these "blank-fields". Fixing L* and 
the normalization to the same value of Grove et al. 2009 
and Hayes et al. 2010 our data is best fit by a substantially 
steeper faint-end slope [a ~ —1.9). As shown by the Cumu- 
lative Luminosity Function (CLF) in Figure [9l this excess of 
Lya candidates is mostly driven by objects with continuum 
detections (blue dotted histogram) , while the CLF for candi- 
dates without continuum detection (red dashed histogram) 
steeply rises only below L < 10*^ erg s~^. The total number 
density above our 90% completeness limit (corresponding to 
Lhya > 4 X 10*^ erg s~^) is n ~ 0.014 cMpc~^, is about a 
factor 5 higher than Grove et al. 2009 and Hayes et al. 2010 
at similar luminosities. 

An enhancement in the number of Lya emitters (LAE), 
i.e. in (p* and possibly in L*, may be expected around a 
quasar simply because this is likely to occupy a high den- 
sity region in the Universe. For instance, Venemans et al. 
(2007) found that LAEs around radio galaxies tend to be of 
a factor 2-5 more numerous with respect to the "field". A 
similar effect may be playing a role in our sample, although 
we note that radio-quiet quasars (like HE0109-3518 ) do 
not generally seem to show a different environment with re- 
spect to "normal", L* galaxies, at least on scales larger than 
0.5 cMpc (e.g.. Groom & Shanks 1999, Serber et al. 2006, 
Falder et al. 2011; see also Hennawi et al. 2006). A very re- 
cent survey by Trainor & Steidel (2012), obtained the same 
result also for hyper-luminous quasars (like HE0109-3518 ) 
at z ~ 2.7, concluding that the very low space density of 
such objects results from an extremely rare event on scales 
<< 1 cMpc and not from being hosted by rare dark matter 
haloes. 

Nevertheless, we also observe a steepening of the LF. 
As previously discussed, fiuorescent boosting should enhance 
the LyQ emission of a gas rich object independent of its SFR 
(or any other combination of parameters that determines the 
Lya luminosity of an isolated galaxy). This constant addi- 
tive boost to the Lya flux will have a larger relative effect 
at low luminosities, thereby steepening the faint end slope 
a. This steepening is the clearest signature of fluorescent 
emission in the luminosity function. 

In order to understand if the observed steepening of the 
LF is compatible with fluorescence, we performed a numer- 
ical experiment with the help of cosmological hydrodynam- 
ical simulations (using the hydro-code RAMSES, Teyssier 
2002) combined with full three-dimensional radiative trans- 
fer for continuum (RADAMESH, Gantalupo & Porciani 
2011) and Lya radiation (Gantalupo et al. 2005; Gantalupo 
et al. 2007). The simulations will be discussed in detail in 
a companion paper (Gantalupo et al., in preparation). In 
Appendix B, we briefiy describe the numerical methods and 
simulation parameters. 

The simulations have been used to produce "mock" ob- 
servations of Lya emitting objects in a 40^ cMpc volume 
of the Universe at 2 = 2.4 around a quasar. The latter has 
been turned on and off in order to simulate the differen- 
tial effect of fluorescent boosting on the LF. The results are 
presented in Fig. 1101 In the "QSO-OFF" simulation we have 
assigned a Lya emissivity to each object proportionally to its 
SFR (see Appendix B) plus the recombination radiation due 
to the UVB background ionizations, i.e. UVB fluorescence 
(several orders of magnitude below our detection threshold) . 
Then, we applied our Lya RT, added the observed noise to 



the simulated image and selected candidates using the same 
SExtractor procedure applied to the real data. The LF ob- 
tained in this way (blue histogram) naturally reproduces the 
observed shape of the LF in "blank-fields" once a normal- 
ization factor niya — 0.13 is applied to take into account the 
fraction of galaxies appearing as Lya emitters and the over- 
density of our simulation box (see Appendix B). We then 
turned on a quasar in the most massive halo in the sim- 
ulation (located at its center) and propagated its ionizing 
photons through hydrogen and helium within the box using 
our radiative transfer code RADAMESH. We then repeated 
exactly the procedure followed for the "QSO-OFF" simula- 
tion, using the same parameters and normalization factor 
of the latter. Thus the only difference between the two sim- 
ulations is the presence of the quasar ionization field. The 
resulting LF in the "QSO-ON" case is shown as a red his- 
togram in Figure [TOl The steepening in the simulated LFs is 
clearly visible and is very similar to the difference between 
the observed LF in the "blank-field" surveys and in our sam- 
ple around the quasar. Within the limitation of numerical 
simulations and observational errors, this provides an ad- 
ditional support for the idea that fiuorescence is playing a 
major role in our sample 

4.3 Luminosity-Distance Relation 

An important prediction for quasar fluorescence is the rela- 
tion for optically thick clouds between the fluorescent Sur- 
face Brightness (SB) and the (three-dimensional) distance 
from the quasar (e.g., Gantalupo et al. 2005). This is be- 
cause the fluorescent boost is proportional to the apparent 
brightness of the quasar as seen by the cloud. Optically thick 
objects whose size is much larger than the ionized region act 
like a "mirror" of the impinging quasar radiation (Gould 
& Weinberg 1996). Their fluorescent SB would thus simply 
scale like the inverse square of their three-dimensional dis- 
tance from the quasar, assuming isotropic and temporally 
constant quasar emission. In reality, however, geometrical 
and kinematical effects can reduce the expected SB by a 
large factor, although this reduction can be parametrized as 
shown in Gantalupo et al. (2005). 

Furthermore, the SB of those gas clouds that have been 
completely ionized by the quasar will be independent of dis- 
tance from the quasar (for as long as they are completely 
ionized). Their luminosity is simply proportional to the in- 
tegral of their hydrogen density squared (see section r5.2p . 
Finally, any local star formation might add a substantial 
Lya flux to the fluorescent component, and this also would 
be expected to be independent of distance to the quasar. 

Despite these caveats, if a significant fraction of opti- 
cally thick clouds ("fiuorescent mirrors") are present, we 
might still expect to find an increase of the average Lya 
SB in our sample closer to the quasar. Unfortunately, we 
cannot obtain directly either the surface brightness or the 
three-dimensional distance from the quasar for most of our 
objects. This is due to their compact nature and the lack, so 
far, of spectroscopic data for the whole sample (a spectro- 
scopic follow-up of the brightest objects will be presented in 
Gantalupo et al., in prep.). However, we can use the total 
Lya luminosity and the projected distance from the quasar 
as proxies for these two quantities, as shown for our sam- 
ple in Fig. 111! Although the statistics are poor, we do find 
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Figure 11. The luminosity of Ly« candidates as a function of 
projected distance from the quasar. The dotted line shows the 
expected fluorescent hya emission produced by an optically thick 
cloud of 1 arcsec'^ area as a function of its three-dimensional dis- 
tance from the quasar. This represents an upper limit for fluo- 
rescent protogalactic clouds without star formation at the quasar 
redshift. The luminosity would scale linearly with object area, but 
the majority of detected objects are compact. One notable excep- 
tion is #1159 (circled), that has an area of ~ 20 arcsec^. The 
dashed line represents the 90% completeness limit of the survey. 
We have used all the objects above this limit (except the very ex- 
tended #1159) to compute the average luminosity per projected 
distance bin (starred open squares). If a significant number of 
optically thick clouds are present and the objects have Lya emis- 
sion boosted by quasar fluorescence, then we expect the average 
luminosity to decrease approximatively with the square of the 
distance from the quasar (assuming no intrinsic relation between 
Lycf emission from star formation and distance from the quasar) . 
Despite the large error bars, there is a marginal detection of the 
expected signal. This is also confirmed by the results of our ra- 
diative transfer simulations presented in figure [T2l 
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Figure 12. Comparison between the average luminosity of ob- 
served hja candidates (black solid line) and the average lumi- 
nosity of the simulated Lya emitters (blue squares for QSO-OFF 
and red dots for QSO-ON simulations) as a function of projected 
distance from the quasar above the same luminosity cut as in fig- 
ure llll Despite large errorbars, the observed trend is very similar 
to the QSO-ON simulation suggesting that the observed signal 
is compatible with fluorescent emission boosted by the quasar 
proximity. 



is turned-on. The trend obtained from the QSO-ON simu- 
lations (red dots) is now in much better agreement with the 
observations (solid line) at all radii. Note that this trend was 
not present without the effect of the quasar radiation (blue 
dots). Within the limitations of our numerical models, this 
is a further support that fluorescent boosting is playing a 
significant role in the Lya emission of the objects detected 
in our survey. 



indeed that the average luminosity above the completeness 
level decreases with the projected distance from the quasar 
(starred open squares). 

As a reference, we have plotted on the same figure the 
expected "mirror" fluorescent Lya luminosity of a optically 
thick cloud with an area of 1 arcsec'^ as a function of three- 
dimensional distance from the quasar (dotted line). This will 
represent the "maximum" emission for a purely fluorescent 
source of this area (e.g. Cantalupo et al. 2005). Although 
the uncertainties related to projection effects and object 
sizes are quite large, it is interesting to note that the ex- 
pected trends are similar to those observed, except in the 
close proximity of the quasar. This is the region where we 
would expect the quasar radiation to be so intense that very 
few clouds would be able to remain optically thick. 

A better analysis that included the size, projection and 
ionization effects can be performed with the help of the ra- 
diative transfer simulations presented in the previous section 
and discussed in detail in Appendix B. In Fig. 1121 we show 
the differential effect of fluorescent boosting on the simu- 
lated projected distance-luminosity relation once the quasar 



5 THE NATURE OF THE SOURCES 

The statistical analyses performed in the previous section 
indicate that the dominant emission mechanism for many of 
our sources is likely to be Lya fluorescence induced by the 
quasar. This mechanism then gives us a unique opportunity 
to directly image, in emission, dense gas independent of any 
associated star formation. This feature distinguishes these 
new observations from previous surveys for Lya emitting 
galaxies or other structures. Moreover, fluorescent boosting 
increases the Lya luminosity of every galaxy, if illuminated 
by the quasar that is gas rich. In principle, this makes in- 
trinsically faint Lya emitters, that would have been below 
our detection threshold, visible in our survey. This raises the 
obvious question: what are the nature and physical proper- 
ties of this new population of sources and how do they differ 
from previously discovered Lya emitters? 

5.1 Proto-galactic clouds and faint galaxies 

As previously discussed, a high Lya EW (e.g., EWo > 240A 
) is the first criterion to exclude objects in which Lya emis- 
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sion is powered (solely) by star-formation. In our sample 
of 98 narrow-band excess objects, we detected 18 sources 
with equivalent width measurements or lower limits above 
EWo > 240A 

Among these 18 objects with demonstrated EWo > 
240 A, aU but six do not have any detectable continuum 
counterpart and therefore quasar fluorescence is likely the 
dominant source of Lya photons. Their Lyo? luminosities, 
which are all around or below 10*^ erg s~^ (see Figure [9l) 
are also fully compatible with fluorescence given their dis- 
tance from the quasar and their compact sizes, as shown in 
Figure 1111 

It should be noted in passing that a further 19 objects 
do not have any continuum detection and could also be en- 
tirely powered by quasar fluorescence. However, the broad- 
band images are not deep enough to require their EWo to 
he above the 240A level. 

It is clear from Figure [9] that the steepening of the cu- 
mulative Luminosity Function with respect to the "blank- 
surveys" involves both objects with and without continuum 
detections - it does not come from the non-continuum de- 
tected objects alone, and is mostly due to the larger num- 
ber of continuum-detected objects that are also presumably 
boosted by quasar fluorescence. The existence of these fluo- 
rescently boosted continuum-detected galaxies is confirmed 
by the EW distributions presented in Figure [7| especially 
when we restrict to the bright sample with Liy^ > 10*^ 
erg s^^. As discussed in section [4.21 our radiative transfer 
simulations naturally reproduce the steepening of the cumu- 
lative luminosity function (see Figure [TCT)) and the variation 
of the average Lya luminosity with distance from the quasar 
(see Figure I12|) once quasar radiation is "turned-on" . It is 
therefore instructive to have a closer look at the physical 
properties of the simulated sources. 

In Figure [T3l we show the relation between the SFR and 
Lya luminosity of the simulated emitters. This has been 
obtained performing SExtractor aperture photometry (af- 
ter adding artificial noise to the image) in the same way as 
with the actual data. In the "QSO-OFF" simulation, the 
only "source" of Lya photons are the stars. The intrinsic 
Lya emissivity is set proportional to the SFR (see appendix 
B for details). In the absence of radiative transfer effects 
the simulated objects in the "QSO-OFF" simulation (blue 
open squares) would thus lie on the relation indicated by the 
black dotted line. It should be noted that our simulations do 
not include absorption by dust particles, and the reduction 
in the recovered luminosity of Lya emitting galaxies is due 
to aperture photometry losses. A significant fraction of Lya 
photons are scattered by the circumgalactic and Intergalac- 
tic medium, forming a diffuse low surface brightness halo. 
This radiative transfer effect was discussed in detail in Can- 
talupo et al. 2005 (see also, e.g., Zheng et al. 2011, Barnes et 
al. 2011). On average, SExtractor aperture photometry only 
recovers about 20% of the total luminosity. This is in agree- 
ment with the recent observational estimate of Steidel et al. 
(2011) from a stacking analysis. Once the quasar radiation 
is turned on, in the simulation, the Lya luminosities change, 
as shown by the solid red dots. Clearly, fluorescent boosting 
makes the simulated galaxies "overluminous" with respect 
to their SFR, increasing the number of detectable objects. 
This is due both to the increased emissivity from recom- 
bination radiation and, to a lesser extent, to the reduced 
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Figure 13. Luminosity versus SFR for the Lya emitters in the 
QSO-OFF (blue squares) and QSO-ON (red circles) simulations. 
In the QSO-OFF simulations, the intrinsic Lya emissivity is set 
proportional to the SFR ( see Appendix B for details). After 
adding simulated noise, the object luminosities are obtained with 
SExactor aperture photometry, similarly to the treatment of the 
actual data. Note that the simulations include the full three- 
dimensional scattering of Lya photons within the circumgalactic 
and intergalactic medium (Cantalupo et al. 2005) but they do 
not include dust absorption. Without radiative transfer effects, 
the Lyo emitters in the QSO-OFF simulation would therefore lie 
on the intrinsic luminosity-SFR relation given by the dotted black 
line, which represents a "maximum" value for the recovered Lya 
emission powered by star formation. Quasar fluorescence in the 
QSO-ON simulation (red dots) boosts many objects above this 
limit, especially at lower SFRs, generating the steepening of the 
CLF observed in Figure [TO] About 20% of sources in the simula- 
tion lying above our nominal detection threshold (dashed line) are 
not associated with star-formation within the simulation (at least 
above our numerical resolution limit of ~ O.OIMq yr~^). They 
are represented by left-pointing arrows at SFR= 0.01. Their rela- 
tive abundances and luminosities are very similar to the observed 
objects without continuum detection in our survey strengthening 
the association between the detected continuum-less sources and 
proto-galactic clouds. 



Lya scattering in the more ionized circum-galactic medium. 
In particular, the simulation suggests that the steepening 
of the observed luminosity function is due to galaxies that 
would otherwise be one order of magnitude fainter in Lya. 
Our survey strategy gives thus the possibility to select and 
study intrinsically faint galaxies at high-redshift at an in- 
comparably deep level. This is similar in some respects to 
the use of gravitational lenses, although in our case only the 
Lya emission is "boosted". In particular, we should have 
reached similar levels to the ultra-deep blind Lya spectro- 
scopic search of Ranch et al. 2008, although a direct com- 
parison is complicated by the very different observational 
techniques. 

An interesting point is that in the "QSO-ON" simu- 
lation, we see sources (the red dots with arrows) that do 
not have any associated star formation at the simulation 
resolution limit for star-formation of about 0.01 Mq yr~^. 
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Figure 14. The relation in tlic simulations between the hya 
luminosity and the maximum HI column density before the quasar 
is turned on (see labels within the panel). The maximum Nhi is 
defined considering all the lines-of-sight that intersect any image 
pixel with detectable Lya emission. Note that we do not simulate 
the formation of H2, and thus the reported HI values will be an 
overestimate, especially at Nhi > lO'^^ cm~^. 



These objects are all below a luminosity of about 10 erg 
s~^ and represent about 20% of the total number of sources 
above our observational threshold. Their luminosities and 
relative abundances are remarkably similar to the observed 
objects that did not have continuum-detection. A closer look 
at these simulated sources reveals that they are compact 
dense clouds, proto-galaxies in the transition phase between 
the IGM and ISM, with peak densities that are not yet large 
enough to currently produce and sustain significant star for- 
mation 0. 

Typically located within dense knots in IGM filaments, 
these clouds will eventually increase their gas density be- 
cause of cosmological gas infall and start to form stars. 
Such a "proto-galactic phase" for gas rich and dense clouds 
might recently have been detected from the study of high 
column density absorption systems such as LLSs (Fuma- 
galli, O'Meara & Prochaska 2011) and among metal-poor 
DLAs (Penprase et al. 2010; Cooke et al. 2011; Cooke, Pet- 
tini & Murphy 2012; but see Carswell et al. 2012). Unfor- 
tunately, the lack of spatial and morphological information 
intrinsic to these one-dimensional absorption studies does 
not allow us to determine whether these clouds are truly 
isolated systems or part of larger galactic gas reservoirs, ei- 
ther circum-galactic (e.g., streams-like features) or within 
galaxies. Given the statistical evidences for fluorescence in 
our sample and the support from simulations, our data pro- 
vides - in emission - strong indications for the existence of 
such compact and isolated proto-galactic or "dark" clouds. 



^ From a numerical point of view, this means that their peak 
density was always below the assumed star formation threshold 
in the simulations of nisM = 1 atoms cm"''. 



opening up a new observational window for the study of the 
early phases of galaxy formation. 

How are these emitters then related to QSO absorption 
line systems? In Figure [14] we show the relation between 
the Lya luminosity of our simulated sources and the maxi- 
mum HI column density before the quasar is turned on. The 
latter is defined as the maximum value among the line-of- 
sight that intersects any pixel with detectable Lya emission. 
We do not simulate the formation of molecular hydrogen, 
therefore the reported HI values are likely to be an overes- 
timate, especially at Nhi > lO'^^ cm~'^. As expected, Lya 
emitting galaxies in the QSO-OFF simulation are only as- 
sociated with the highest Nhi systems, i.e, Nhi > 5 x 10^° 
cm~^, above the DLA threshold (2 x 10^° cm"^). This is in a 
sense by construction, since we only allowed star-formation 
to take place above a given density threshold (1 atom cm""') 
following the Schmidt-Kennicutt law (see Appendix B). In- 
deed, this density threshold translates into a critical column 
density (given by the Toomre instability criterion ) that in 
local galaxies is observed to range between 5 x 10^" cm""^ 
and 2 x 10^^ cm"^ (Kennicutt 1998). 

Once the quasar is turned on, we notice that the overall 
boost in luminosity still seems to correlate with the initial 
HI column density of the quasar-off simulation. This is be- 
cause higher column density systems are typically larger, 
more massive and denser. The proto-galactic clouds with- 
out star-formation start to appear at values of Nhi around 
or below the simulation critical SF threshold (5x 10'^'' cm""^). 
In particular, the simulation suggests that the objects fluo- 
rescently illuminated by the quasar have typical Nhi - when 
the quasar is off - one order of magnitude lower than the 
sources detectable by similarly deep "blank-field" surveys. 
Given our current sensitivity limits, we are able to detect 
only the high-end tail of the Nhi distribution associated 
with these clouds. Deeper surveys, e.g. with integral field 
spectrographs, could reveal many more such systems. 



5.2 The gas mass and star formation efflciency of 
protogalactic haloes 

If we can assume that there is negligible Lya emission pro- 
duced by internal star formation, then we can directly con- 
vert the Lya luminosity produced by HII recombination into 
an estimate of the total HII mass. Indeed, in the absence 
of dust absorption and HI collisional excitations (in other 
words, if the gas is mostly ionized), then the total Lya lu- 
minosity is simply given by (e.g., Cantalupo et al. 2008): 

Liyc, = / rio • nnii • aiy^{T)huiyadV, (2) 
Jv 

where Uc, nmi, a{^^^(T) and huiya are, respectively, the elec- 
tron number density, the proton number density, the effec- 
tive Lya emission coefficient (see Cantalupo et al. 2008), 
and the Lya photon energy. Assuming a uniform tempera- 
ture and ionization fraction within the cloud, we can rewrite 
eq[5]as: 

Liya = XcOiya{T)hiyiya / nliiidV, (3) 



© RAS. MNRAS 000.[Tlf26l 



18 



10" 



Simulation QSO-ON 

□ SFR>0.1 Mg/yr 
fo SFR<0.1 Mg/yr 



no"' 




10"' 



Figure 15. The relation between the cloud gas mass and hya 
luminosity within our simulation with the quasar turned on. We 
have included only gas with T < 5 X 10*K. The solid line is 
the analytical expected relation (eq[8} for fully ionized objects 
without star formation, derived from the halo cosmological scal- 
ing laws and assuming a fixed gas temperature (see text for de- 
tails). This relation gives a good approximation for the expected 
mass-luminosity function of fluorescently illuminated gas clouds 
with absent or little star-formation, as derived from our radiative 
transfer simulation. Both analytical arguments and the numeri- 
cal simulation suggest therefore that the gas mass of the detected 
protogalactic candidates would, if it is mostly ionized, be of the 
order of 10^ Mq or larger. 



where, Xo is the electron fraction 1.16 if H and He are 
fully ionized). Using the definition of clumping factor: 



C ■- 

we thus obtain: 



< nun > 



[Jy nniidV] 



(4) 



(5) 



For a spherical, fully ionized cloud with radius R, gas mass 
Mgas and fixed temperature in the range lO'^K < T < 
W'^'^K, eqO becomes: 



1.4xlO**'^erg/s-C- 



2 • lO^K 
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n -3 r 



pkpc 



109 M© 



(6) 

where pkpc denotes physical kpc. With the knowledge of the 
cloud size we could therefore derive its gas mass from the 
Lya luminosity. 

Unfortunately, our proto-galactic candidates are essen- 
tially unresolved in our ground-based image. In absence of 
this information, we can assume that the gas is settled into a 
disk at the center of its dark matter halo and derive a proxy 
for the object size from the cosmological scaling relations be- 
tween virial radius and mass (e.g.. Mo, Mao & White 1998). 
In particular, assuming that the disk size and gas mass are 
fixed fractions (respectively, fr and fg) of the dark-matter 
halo virial radius and total mass, we have (at z = 2.4): 



R ~ 2.7pkpc 
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(7) 



where we have used the same values of the disk to total 
mass ratio (fg) and disk to virial radius ratio (fr) as in 
Mo et al. (1998). It should be noted for consistency that 
eq[7] implies that a disk with Mgas ~ 10^ AIq would have a 
size of less than 0.7" at this redshift and therefore it would 
appear unresolved in our image. Using eq[7] as a proxy for 
the cloud size (and dropping the dependency on fg and fr 
for simplicity) we can thus rewrite eq[6] as: 



1.4 X W'Pl 







lO^'^erg s~i 



T 



2 X W^K 



■C" 



(8) 



We stress that the mass derived with this approach is rep- 
resentative of the cloud gas mass only for systems that are 
highly ionized by the quasar radiation and with a given tem- 
perature. In the other cases, eq[8] underestimates the value 
of Afgas, unless there is a significant Lya emission contribu- 
tion from HI coUisional excitation, e.g. from cooling. How- 
ever, highly ionized, optically thick objects may be distin- 
guished from cooling clouds using the Lya spectral and SB 
profile (see, e.g. Cantalupo et al. 2005, Dijkstra et al. 2006, 
Cantalupo et al. 2007) or with the help of other diagnostic 
lines such as Hell 1640A (e.g., Yang et al. 2006, Scarlata 
et al. 2009) . Future spectroscopic follow-up and space-based 
imaging of our sources may therefore help to understand if 
eq[8]is applicabile in these objects. In the current absence of 
this information, however, we can use our radiative transfer 
simulation to have an idea of the expected mass-luminosity 
relation for a variety of objects with different ionization frac- 
tions and temperature distributions. 

In Fig. [TS] we show the relation between the cloud gas 
mass and its Lya luminosity in our radiative transfer sim- 
ulations, once the quasar is turned on. We have considered 
only gas with T < 5 x 10^ K - the main component produc- 
ing Lya photons from recombination - and lying within the 
aperture radius given by SExtractor photometry. We divide 
the sample according to the SFR of the object. For clouds 
with SFR< O.IMq yr~^ (red circles), the Lya luminosity is 
mostly powered by quasar radiation. These systems appear 
to follow well the analytical trend with mass. As expected, 
eq[8] tends to underestimate the total gas mass, although we 
note that a higher temperature than 2 x 10* K would help 
to recover the value obtained in the simulations. 

From eq[8] and with the help of the radiative trans- 
fer simulation, we can therefore derive an estimate of the 
gas mass of those Lya sources without continuum detec- 
tion and with EWo > 240A. Their typical Lya luminosity 
are in the range 0.5 — 1 x 10*'^ erg s~^ implying therefore 
a gas mass of around Mgas ~ Mq, or higher if they are 
not mostly ionized by the quasar radiation. We can derive 
a constraint on the average SFR of these systems from the 
non-detection in the V-band stack (see FigU} at the level 
of V{la) ~ 30.3. We convert this observed magnitude limit 
(corresponding to rest-frame wavelength of about I6OOA) 
into a SFR upper limit using Starburst99 (Leitherer et al. 
1999) and Oti-Floranes & Mas- Hesse (2010). Assuming an 
extended burst of duration 250Myr, Salpeter IMF (1-100 
Mq) and E{B - V) = we obtain a SFR< W'^Mq yr"^ 
This implies a Star Formation Efficiency (SEE) for these 
objects of SFE= SFR/Mgas < 10"" yr'V 

This very low SEE is several times lower than the typ- 
ical SFE of dwarf galaxies in the local Universe (e.g., Geha 
et al. 2006) . It is also 200 times lower than the SFE of typ- 
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Figure 16. Continuum subtracted NB image of #1159 reveal- 
ing the extended, filamentary emission surrounding the central 
source. Contours show isophotal surfaces (Itr intervals, starting 
from 2(7). Size and luminosity of this emission are compatible 
with the circumgalactic medium illuminated by the quasar (see 
text for details). Note that 8" corresponds to about 66 physical 
kpc at z=2.4, i.e. the virial diameter of a halo with total mass 
A/ ~ 4 X 10"Mq. 



ical main sequence galaxies at 2 ~ 2 which has been es- 
timated to be around SFE~ 2 x 10~^ yr~^ (Daddi et al. 
2010, Genzel et al. 2010). Although we currently have only 
upper limits on the cloud sizes, it is interesting to estimate 
where our sources would lay with respect to the Kennicutt- 
Schmidt (KS) relation (Kennicutt 1998) between SFR per 
unit area (Ssfr) and gas surface density (Egas). From the es- 
timated gas mass Mgas ~ 10® M©, the SFR upper limit and 
assuming that our objects follow the cosmological scaling 
relation between mass and radius (eq. [7]), we would obtain 
an average Egas ~ 44 Mq pc~^ and Esfr < 4 x lO"* Mq 
yr~^ Kpc"'^. This value of Esfr is about 100 times lower 
than the expected value from the local KS relation but 
it is consistent with the metallicity-dependent models of 
Krumholz, McKee and Tumlinson (2009) and Gnedin & 
Kravtsov (2011) that predict a substantial steepening of the 
KS relation at Egas ~ 50 — IOOA/0 for very low metallicity 
gas {Z < O.OIZ0) at high redshift (see also Rafelski, Wolfe 
& Chen 2011). 

5.3 Cold Circum-Galactic Medium 

A different category of objects is represented by the more 
luminous, continuum-detected sources with EW> 240A, 
which show evidence of extended emission (left-hand pan- 
els in Figure [5]). These are among the brightest Lya sources 
in our sample and therefore it is not surprising that we de- 
tected some stellar continuum associated with them. Their 
extended emission has a surface brightness at a level of 
0.5 — 1 X 10~^^ erg s~^ cm~^ arcsec"^, a value compati- 
ble with the quasar induced fluorescence suggested by the 
high EW. Their luminosities and low number density sug- 
gests, on the basis of our simulations, an association with 
the most massive systems in our surveyed volume, i.e. with 
objects with halo masses between lO" and 10^^ Mq. Sev- 
eral theoretical and numerical studies have suggested that 



these haloes are able to retain dense and cold gas reservoirs 
in the form of filaments in the circum-galactic medium (e.g., 
Dekel & Birnboim 2006, Ocvirk & Teyssier 2008, Dekel et 
al. 2009, Keres et al. 2009, Fumagalh et al. 2011, van de 
Voort & Schaye 2012). At the densities and temperatures 
predicted by these models, these circum-galactic filaments 
could be detectable in fluorescent emission in the proximity 
of a bright quasar. Note that this is a different mechanism 
with respect to Lya cooling emission - suggested as a possi- 
ble origin for extended Lya structures like Lya blobs, away 
from bright quasars (e.g., Haiman, Spaans & Quataert 2000, 
Fardal et al. 2001, Dijkstra & Loeb 2009, Goerdt et al. 2010; 
but see Furlanetto et al. 2005, Faucher-Giguere et al. 2010; 
Rosdahl & Blaizot 2011). Mostly due to HII recombinations, 
quasar fluorescent emission is much less sensitive to the gas 
temperature state (see e.g.. Fig. 1 in Cantalupo et al. 2008) 
and therefore it is much easier to estimate the expected Lya 
surface brightness (see Appendix B for further discussion). It 
is intriguing that one of the brightest objects in our sample, 
#1159, indeed shows some evidence for possibly filamentary 
emission features (Fig. I16p extending over 8" in the sky. At 
z — 2.4, this corresponds to a projected size of 66 physi- 
cal kpc, i.e. the virial diameter of a halo with total mass 
M ~ 4 X 10^^ M0. Without kinematical information, it is 
however difficult to determine if these extended features are 
associated with infall or outflows. 

Other mechanisms can produce fllamentary structures. 
In fact, in our simulations, the most common sign of fll- 
amentary CGM emission is associated with tidal features, 
e.g., dense gas stripped out of satellite galaxies, rather than 
the cold infall of gas. This could be the origin of the extended 
emission observed in our images. In particular the observed 
"bridge" between two of our bright extended sources, #2411 
and #2433 might fall into this category (see also Ranch et 
al. 2011). Further theoretical and observational studies fo- 
cusing on these sources, e.g., with the help of Integral Field 
Spectroscopy, might help reveal the detailed properties of 
the detected extended emission. 

While much remains to be learnt about these systems, 
the key point is that our narrow-band survey has demon- 
strated the potential of flnding and mapping the extended, 
cold circum-galactic medium around galaxies thanks to flu- 
orescent boosting of the emission. This mechanism opens up 
a new window for the study of this important gas reservoir of 
high-redshift galaxies and to unraveling how galaxies receive 
their gas. 

6 UNCERTAINTIES AND LIMITATIONS 

6.1 Foreground and AGN contamination 

What fraction of our Lya candidates could be foreground 
line emitting galaxies and how do this affects our results? 
Typically, the primary sources of contamination for high- 
redshift NB Lya surveys are foreground galaxies with strong 
[Oil] 3726-3729A emission line doublet. For Lya surveys at 
2 < 2 < 3, corresponding to 0.1 < 2 < 0.3 for [Oil], the 
typical fraction of foreground contamination is typically less 
than 5% (e.g., Venemans et al. 2007), because of the much 
smaller volume and smaller EW associated with [Oil] emit- 
ters. From our data, we cannot directly estimate the fore- 
ground contamination in our sample because of the lack of 
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spectroscopic confirmation for the wliole sample. The spec- 
troscopic follow-up of about twenty among the brightest ob- 
jects in our sample (Cantalupo et al., in prep.) clearly shows 
that they are Lya emitters but the intrinsic faintness for the 
rest of our sample makes their direct confirmation more dif- 
ficult. 

Although foreground contamination could increase the 
number of detectable objects (affecting the luminosity func- 
tion), it cannot produce our skewed EW distribution to- 
wards higher values. On the contrary, given their smaller 
EW, foreground objects would be expected to produce the 
opposite effect. Moreover, they would not show any correla- 
tion between their luminosity and the position of the quasar. 
Therefore, we believe that the presence of foregrounds in our 
sample does not affect our conclusions. 

Unlike star-powered galaxies, active galactic nuclei 
(AGN) are able to produce EW well above the canonical 
240A value and might thus be "confused" with fiuores- 
cent emission. AGN can be identified by spectroscopic fea- 
tures (e.g., broad emission lines and the presence of the CIV 
doublet) or from X-ray coincident detections. Previous sur- 
veys have shown that AGN contamination is always very 
low, of the order of few percent (e.g., Guaita et al. 2010). 
It is in principle conceivable that there might be an in- 
creased number of AGN around our ultra-luminous quasar. 
Although clustering effects would increase the number of 
galaxies around the quasar, to our knowledge there are no 
indications in the literature for the increase, by a factor of at 
least five, in the AGN fraction in such environments, which 
would be necessary to explain our results in terms of AGN. 
This is corroborated by our preliminary spectroscopic con- 
firmation for the bright part of our sample (Cantalupo et 
al., in prep.): only one object out of about twenty shows the 
presence of CIV emission. Therefore, although AGN might 
contribute to increasing the number of high EW objects in 
our sample, it is very unlikely that they could account for 
observed difference with the "blank-fields" . 

6.2 Collisionally excited Lya 

Lya produced by coUisional excitation via energetic elec- 
trons, commonly called "Lya cooling" , is a mechanism that 
is in principle able to produce strong Lya emission with- 
out associated star formation, i.e. similarly to quasar flu- 
orescence. Unfortunately, collisionally excited Lya emission 
has an exponential dependence on gas temperature (see e.g., 
Fig.l in Cantalupo et al. 2008), which makes accurate pre- 
dictions within current hydro-simulations very challenging 
(see Appendix B for discussion). 

For the primary purpose of this study - finding proto- 
galactic clouds and direct imaging of the cold CGM - it 
does not really matter if Lya photons are produced by 
quasar fiuorescence or by cooling due to gravitational col- 
lapse. From an observational point of view, Lya cooling 
should be equally detectable away from bright quasars. If 
our sources are in fact powered by cooling instead of fluo- 
rescence, then the "blank-field" surveys at 2 < 2 < 3 should 
already have detected such objects, and we should not see 
any particular difference between, e.g. our EW distribution 
and the one obtained by these surveys. On the other hand, 
cooling Lya emission could explain why "blank-field" sur- 
veys do detect some objects with high EW away from bright 



quasars, including the large Lya blobs. These sources would 
also then be present within our survey. However, given their 
low space density (e.g., Matsuda et al. 2004) we would ex- 
pect that they are responsible for only a small fraction of 
our sample of high EW sources. 

6.3 Quasar eflfects 

The ionizing radiation from the quasar has the positive ef- 
fect, for us, of increasing the ionization state, and thus the 
recombinant Lya emission, of the gas. However, the same 
radiation could also have an important effect on the galaxy- 
formation process itself around the quasar, e.g., reducing the 
gas cooling rate (see Cantalupo 2010) and thus effectively 
limiting gas accretion onto (proto-)galaxies and suppressing 
their SFRs. In principle, from the continuum properties of 
the sources in our survey it would be possible to search for 
a general star-forming deficit relative to the "blank-fields". 
However, the intrinsic faintness of our objects plus uncer- 
tain clustering effects would make this task very difficult in 
practice. 

Nevertheless, let us assume for the moment that the 
galaxy-formation process would be "slowed-down" for a gas 
cloud while it is "illuminated" by the quasar, e.g. if the 
gas is completely ionized. In this case, the overall effect on 
the galaxies would depend on the quasar emission history 
compared to the galaxy-formation timescales. Within the 
current paradigm of quasar black-hole accretion the quasar 
luminous phase is expected to have a relatively short du- 
ration (e.g. a few tens of Myrs) followed by a much longer 
phase (e.g., a few Gyrs) of lower ultraviolet output (e.g., Di 
Matteo et al. 2008; see also Martini 2004 for a review). The 
dense gas, which should be associated with quasar fluores- 
cent emission has a recombination time scale is very much 
shorter, of the order of 10^ — 10® yrs, than the quasar life- 
time, and of cosmological timescales. This means that any 
previous bright phase of the quasar will be quickly "forgot- 
ten" once the quasar turns off, and proto-galactic clouds 
and galaxies should be able to recover their "normal" evo- 
lution before the current phase. At the same time, the star 
formation timescale probed by our deep V band image (cen- 
tered on A ~ I6OOA) is at least 200 Myrs (e.g., Leitherer et 
al. 1999), one order of magnitude larger than the expected 
quasar luminous phase. Therefore it is unlikely that the cur- 
rent bright phase is the cause of the apparent lack of star 
formation in our proto-galactic clouds candidates. Of course, 
we cannot exclude that our quasar had a very peculiar emis- 
sion history. Future observational campaigns will be crucial 
to increase the the sample to a larger number of quasars and 
therefore reduce these uncertainties. 



7 SUMMARY AND CONCLUSIONS 

We have presented a deep survey for fluorescent Lya emis- 
sion in a large cosmological volume (~ 5 x 10"^ comov- 
ing Mpc) centered on the hyperluminous 2 = 2.4 quasar 
HE0109-3518. With the help of a custom-made, narrow- 
band filter mounted on the VLT/FORS spectrograph, we 
were able to photometrically select 98 Lya candidates down 
to a detection limit of about 4 x 10~^* erg s~^ cm~^ (cor- 
responding to Lhya ^ 2 X 10^^ erg aX z = 2.4). 
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The properties of these Lya sources in terms of their: 
i) equivalent width distribution, ii) shape of the luminosity 
function and, iii) the variation of their luminosity with pro- 
jected distance from the quasar, all indicate that our sample 
is consistent with having much of the Lya emission originat- 
ing in fluorescent reprocessing of quasar radiation. 

In particular, the main evidence supporting the fluores- 
cent origin of our sample are: 

• The large number of objects with rest-frame EWo > 
240 A (18 sources out of 98), commonly assumed to be the 
limit for Lya emission powered by Population II stars (e.g., 
Schaerer 2002). Independent of luminosity, the fraction of 
such objects appears to be about 20%, which is about a 
factor of four larger than "blank-field" surveys (e.g., Ciar- 
dullo et al. 2012, Grove et al. 2009) or surveys around radio- 
galaxies (Venemans et al. 2007) at similar redshifts. More- 
over, the whole shape of our equivalent width distribution is 
skewed toward higher values, with an exponential scale that 
is about double that found in the "blank-fields" . 

• The steeper faint-end of our LF with respect to the 
"blank-fields" surveys with a corresponding increase in the 
source number density by a factor of about five compared to 
these studies. The steepening of the LF is well reproduced 
by a radiative transfer simulation of fluorescent emission 
around quasars. The simulation suggests that the difference 
with the "blank-field" surveys may be fully explained by 
fluorescent boosting of gas-rich but intrinsically faint Lya 
emitters. 

• The increase in the average luminosity of the sources 
at smaller (projected) distances from the quasar. Although 
this effect could also be caused by clustering of bright sources 
around the quasar, our simulations suggest that the observed 
trend is very similar to the theoretical expectations for flu- 
orescent emission. 

Fluorescent emission gives us a unique opportunity to 
directly image, in emission, dense, high-redshift gas inde- 
pendent of any associated star formation. In particular, we 
identify three categories of objects associated with fluores- 
cent emission in our study: 

• Proto-galactic or "dark" clouds, i.e. gas-rich objects 
with little or no associated star-formation. Candidates for 
these sources are represented by the 12 objects without de- 
tectable continuum and with lower limits on their EW higher 
above 240A . They are typically compact or unresolved with 
L < 10*^ erg s~^. A V-band image stack of these objects did 
not show any continuum detection above la {V ~ 30.3), im- 
plying a combined constraint on their EW of EWo > 800A 
(la). This effectively rules out Lya powered by internal star 
formation. A comparison with our radiative transfer simu- 
lations suggests that these sources would have maximal HI 
column densities around or slightly below the DLA thresh- 
old (2 X 10^" cm~2) if they were observed away from the 
quasar. These may therefore be the analogues of absorp- 
tion systems which have also been suggested to probe such 
a "proto-galactic phase" (e.g., Fumagalli et al. 2011, Cooke 
et al. 2011). Observing these systems in emission - rather 
than in absorption - allows us to readily discriminate be- 
tween truly isolated clouds and sub-components of larger, 
galactic or circum-galactic gas reservoirs. Our narrow-band 
imaging provides indications for the existence of compact 



and isolated proto-galactic or "dark" clouds opening up a 
new observational window for the study of the early phases 
of galaxy formation. Assuming that these objects are al- 
most fully ionized by the intense quasar radiation, we can 
directly derive an estimate or their total gas mass, namely 
Afgas ~ 10® Mq. From the non-detection in the V-band 
stack, we can furthermore estimate an upper limit on their 
average star formation rate of SFR< 10~^Mq yr~^. This 
implies that their star formation efficiency is lower than 
SFE< 10"^^ yr~^ (or, equivalently, their gas consumption 
time is larger than 100 Gyr), several times below the SFE 
of gas- rich dwarf galaxies in the local Universe (e.g., Geha 
et al. 2006) and two hundred times lower than in typical 
massive star-forming galaxies at these redshifts (Daddi et 
al, 2010, Genzel et al 2010). These objects are therefore the 
best candidates for proto-galactic clouds or "dark" galaxies 
at high-redshift suggested by recent theoretical studies (e.g., 
Bouche et al. 2010, Gnedin & Kravtsov 2010, Krumholz & 
Dekel 2011, Kuhlen et al. 2012). 

• Intrinsically faint Lya emitting galaxies fluorescently 
boosted above our detection threshold. As suggested by our 
radiative transfer simulations, these sources are the main 
component producing the observed e-folding scale increase 
in the EW distribution and the steepening of the LF. In 
particular, their implied SFRs are about one order of mag- 
nitude smaller than the values we could have detected in 
absence of quasar fluorescence. 

• Extended, cold circum-galactic medium emission. At 
least three of the six bright sources with EWo > 240A and 
continuum detection in our sample show extended, possi- 
bly filamentary fluorescent emission. In particular, one ob- 
ject (#1159) presents signs of several filaments connected to 
a bright central source. The size and morphology of these 
structures are compatible with the expectations from theo- 
retical and numerical models for cold streams feeding high- 
redshift galaxies (e.g., Dekel et al. 2009). Our current data 
and the lack of kinematical information cannot however rule 
out a different origin for this cold and dense circum-galactic 
gas, e.g., tidal stripping from gas-rich satellite galaxies (as 
commonly seen in our simulations). 

After several years of attempts to detect fluorescent Lya 
emission, our results demonstrate the potential of a carefully 
designed (very-)narrow-band survey around bright quasars 
to discover and study this new category of sources. Sev- 
eral important factors have contributed to this success. As 
we have shown, a large survey volume, comparably deep 
narrow-band and continuum imaging, as well as sophisti- 
cated numerical models including radiative transfer effects, 
are fundamental to detect, recognize and interpret fluores- 
cent emission. Still, the success of such a fluorescence sur- 
vey might depend on the emitting history and properties of 
the source that provides the "illumination", i.e. the quasar. 
In this respect, our survey suggests that for HE0109-3518 
at least, the proximity effect extends on much larger an- 
gular scale than just our line-of-sight, the intrinsic limit 
of absorption-line studies. Object by object variations, in- 
cluding clustering and cosmic variance, might however play 
an important role and observational campaigns targeting a 
larger number of quasars will thus be important in the fu- 
ture. 

By stimulating fluorescent emission in neutral gas, the 
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quasars are literally shedding light on the hitherto obscure 
early stages of galaxy-formation and on the way in which 
galaxies acquire their gas. With the present study, we make 
a crucial step towards revealing and understanding how to 
take full advantage of this unique opportunity. 
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APPENDIX A: LIST OF CONTINUUM DETECTED OBJECTS WITH EWo < 240A 



Table Al: Position and properties Lya candidates with continuum detection and EWo < 240A (see Table 2) 
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APPENDIX B: RADIATIVE TRANSFER SIMULATIONS 

In this section, we briefly describe the hydrodynamical and Radiative Transfer (RT) simulations used throughout the paper. 
A detailed description will be provided in a separate paper (Cantalupo & Haehnelt, in preparation). 

We have performed cosmological Adaptive Mesh Refinement (AMR), hydrodynamical simulations using the publicly 
available code Ramses, version 3.03 (Teyssier 2002). The computational domain is composed by a series of four nested boxes 
at increasing, root-grid resolution. The largest box has a linear size of 40 comoving Mpc, corresponding roughly to the FWHM 
of the NB filter used in our FORS2 observation at z = 2.4, while the high-resolution box has a linear size of 10 comoving 
Mpc, corresponding roughly to the FOV of our FORS2 observation at 2 = 2.4. The equivalent base-grid resolution in the 
high-resolution region corresponds to a (1024'') grid with a dark-matter particle mass of about 1.8 x 10® M©. We allowed 6 
additional refinements level, reaching a maximum spatial resolution of about 0.6 comoving kpc, i.e. about 180 proper pc at 
z = 2.4. The box has been centred on the most massive halo at z = 2.4, Mdm ~ 10^"^ Mq, as a representative host of a 
luminous quasar. We include photo-ionization from the cosmic UV background (Haardt & Madau 2012) using a self-shielding 
above a critical density as derived by Schaye (2001). In each cell, gas is converted into star particles following a Schmidt law 
with a density threshold of 1 atom cm"^ and a star formation timescale of 3 Gyr (see, e.g. Dubois & Teyssier 2008). We 
include metal cooling, metal enrichment and supernova feedback as described in Dubois & Teyssier 2008. 

To produce the "quasar-on" models, the results of the hydrodynamical simulation have been post-processed with our RT 
code RADAMESH (Cantalupo & Porciani 2011) to include the time-dependent ionization effects due a bright quasar with 
Lyman Limit luminosity of Ln = 5 x 10"'^ erg s~^ and spectral slope of a = —1.7 (Teller et al. 2002). We include hydrogen 
and helium ionization and we evolved the temperature of the gas consistently with the ionization state using a large number 
of frequency bins (50) for the ionizing radiation. We assumed a quasar age of 20Myr, large enough to illuminate the whole 
simulation box. 

In order to obtain simulated Lya images we post-processed the "quasar-off" and "quasar-on" runs with our three- 
dimensional Lya Monte Carlo RT (Cantalupo et al. 2005) including the gas velocity field and inhomogeneous temperature 
distribution. As a source of Lya radiation we included HII recombination and internal star formation, while we have decided 
to exclude contribution from HI coUisional excitations. This is due to the fact that a proper modeling of this emission for 
dense systems would have required the coupling of hydro and RT, given its strong temperature dependence for a photo-ionized 
gas (e.g., Cantalupo et al. 2008; see also, e.g., Furlanetto et al. 2005, Faucher-Giguere et al. 2010 and Rosdahl & Blaizot 2011 
for further discussion), currently beyond the capabilities of the large majority of numerical studies. We note however that 
the contribution from HI coUisional excitations is likely negligible if the gas is highly ionized, as expected in proximity of 
an ultraluminous quasar. Partially ionized systems, e.g., optically thick clouds to the quasar radiation could present instead 
significant contribution from HI coUisional excitations. For these sources, the simulated Lya emission, in absence of dust, 
needs to be considered a conservative estimate. 

We computed the Lya emissivity from the simulated Star Formation Rate (SFR) converting the Ha-SFR relation of 
Kennicutt (1998) to a Lya-SFR relation using the Case B line ratio (e.g., Osterbrock 1989), obtaining: SFR= 9.1 x 10~''3Liya 
erg s~^. We do not include dust absorption, likely important only for the most massive Lya emitters powered by internal star- 
formation, while the "typical" Lya emitter seems to contain very little dust (e.g., Guaita et 2011). Nevertheless, continuum- 
selected samples of star forming galaxies at z 2.5 (e.g., Steidel et al.2011) show that only about half of the galaxies 
present Lya in emission and only 20% of them have EWo > 20A, the standard selection cut. In other words, not every 
star-forming galaxy is a Lya emitter and we must take this into account when deriving the simulated Luminosity Functions. 
In particular, when we computed the simulated Luminosity Function in the "quasar-off" model (where Lya is only powered by 
star formation), we have used a normalization factor niy^ = 0.5/3.7 = 0.13, where the factor 3.7 is necessary to renormalize for 
the overdensity of the simulated high-resolution region. Note that, given the large uncertainties related to the Lya fraction and 
its evolution with galaxy luminosity (e.g.. Stark et al. 2010), niya is essentially a free parameter in our model. It is remarkably, 
however, that with the adopted value we are able to obtain a simulated LF very similar to the observed LF by Grove et al. 
(2009) and Hayes et al. (2010). In the "quasar-on" simulation we use the same rescaling factor, since we are interested in the 
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change of slope rather than normahzation. Also in this case, however, the adopted value of n\ya - combined with the steepening 
of the LF - seems to produce the right amount of sources, regardless of the possible difference in overdensity between quasar 
and random fields. Note that in the "quasar-on" simulation we have assumed isotropic quasar emission, while in reality only 
a fraction of the volume may be illuminated by the quasar. This parameter would be degenerate with field overdensity. In 
other words, the adopted value of niyo, for the fluorescent sources may also be consistent with a cosmic volume 4 times more 
dense than a random field, if the quasar is only "illuminating" one fourth of such volume. 
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